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This  report  examines  point  detects  in  ice:  molecular  d^ects  (vacancies  and 
interstitials),  protonic  point  delects  (ions  and  Bjemjm  defects)  and  electronic 
point  defects  (sotvaled  electrons  and  radicals).  Experimental  results  and 
theoretical  models  ofthed^ts' atomic  and  electronic  structures,  energies  and 
mobilities  are  reviewed.  Special  attention  is  given  to  the  results  of  studies  from 
the  last  two  decades.  Among  the  experimental  techniques  under  consideration 
are  nuclear  magnetic  resonance.  X-rays  topography,  electron  spin  resonance, 
measurements  of  Ice  conductivity  and  dielectric  permlttivily,  spectra  of  optical 
absorption  and  photoiuminescence,  and  measurements  of  diffusion  co^cients 
of  the  oxygen  and  hydrogen  isotopes  in  ice.  Ab  initio  calculations,  computer 
simulations,  Jaccard's  theory  and  solilon  models  are  used  to  discuss  the 
defects'  structures  and  their  properties. 
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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  euround  the  world  deal 
with  ice,  snow  and  permafrost.  They  are  scientists,  educators,  engineers,  navigators, 
meteorologists  and  others.  While  a  small  fraction  of  these  people  c<»ttribute  to  the 
knowledge  base  in  ice  physics,  all  of  them  use  knowledge  from  it  frequently. 
Moreover,  successful  applied  research  is  based  upon  fundamental  science — one 
more  reason  for  ice  specialists  to  have  a  textbook  on  ice  physics  on  their  desks. 

The  first  modem  ice  physics  text  was  Fletcher's  book  on  The  Chemical  Physics  of  Ice 
(1970).  Fletcher's  book  is  in  typical  textbook  format:  it  is  reasmably  brief  and  easy 
to  understand.  He  touched  on  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in 
1974.  Hobbs  considered  almost  <ill  of  the  basic  aspects  of  ice  as  understood  at  that 
time.  Moreover,  he  described  and  compared  several  (sometimes  opposing)  view¬ 
points.  This  fundamental  and  rather  large  (837  pages)  book  is  commr^y  known  as 
the  "Ice  Bible"  by  specialists  in  the  field.  In  1974  and  1975,  two  CRREL  Monographs 
on  ice  were  produced  by  John  Glen.  These  were  briefly  and  clearly  written  and 
reviewed  almost  all  ice-related  subjects.  This  work  was  (and  in  some  respects  still  is) 
a  magnificent  introduction  to  ice. 

Finally,  in  1981  Maeno  wrote  a  simple  popular  book  for  the  express  purpose  of 
attracting  people's  attention  to  the  subject. 

During  the  past  20  years,  a  significant  amount  of  new  experimental  and  theoreti¬ 
cal  work  has  appecired,  dramatically  changing  our  views  on  ice  physics.  As  a  result, 
we  are  now  able  to  formulate  physical  laws  using  more  simple  and  direct  methods. 
We  have  found  some  of  the  physical  modeb  used  in  the  past  to  be  completely  wrong. 
The  physics  of  ice  is  a  much  Iwtter  developed  subject  than  it  was  20  years  ago. 

For  the  above  reasons,  we  feel  the  time  is  ripe  for  a  contemporary  book  on  ice 
physics,  incorporating  the  known  and  proven  with  almost  20  years  worth  of 
material  not  covered  by  previous  works. 

We  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader 
to  be  a  uniquely  educated  person.  It  is  our  intent  to  present  the  material  in  such  a  way 
that  ciny  reader  attracted  by  the  title  Ice  Physics  wiU  be  able  to  comprehend  it.  This 
is  quite  difficult  for  a  book  dedicated,  not  to  a  particular  field  of  knowledge,  but  to 
a  specific  material.  Indeed,  for  ice  it  means  we  have  to  consider  a  wide  variety  of 
subjects,  including  quantum  chemistry,  solid  state  physics,  the  theory  of  elasticity, 
ionic  conductivity,  synchrotron  x-ray  topography,  crystal  growth,  fte  physics  of 
surfaces  and  more. 

The  primcuy  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing 
scientific  accuracy.  Experimental  facts,  physical  ideas  and  theories  will  be  strongly 
organized  and  bound  together  cohesively.  The  reader  will  be  introduced  to  a  wide 
variety  of  material  on  a  step  by  step  basis.  Then  the  picture  will  be  whole. 

To  accelerate  materuils  publication,  this  book  will  appear  first  in  the  form  of  a 
series  of  joint  CRREL-Dartmouth  reports,  later  to  be  published  in  CRREL's  Mono¬ 
graph  series,  on: 

1.  The  structure  of  ordinary  ice 

Part  I:  "Ideal"  structure  of  ice.  Ice  crystal  lattice 
Part  11:  Defects  in  Ice 

Volume  1:  Point  defects 

Volume  2:  Dislocations  and  planar  defects 

2.  Electrical  properties  of  ice 

Part  I:  Conductivity  and  dielectric  permittivity  of  ice 
Part  H:  Advanced  topics  and  new  physical  phenomena 

3.  Optical  properties 

4.  Electro-optical  effects  in  ice 


iii 


5.  Themial  properties 

6.  Mechanical  pn^wrties  ice.  Elasticity  and  anelastk  lelaxatkm.  nastk 
properties.  Fracture  ice 

7.  Etectromechankal  effects  in  ice 

8.  Surface  of  ke 

9.  Other  forms  of  ice  and  their  properties 

10.  keinsp^e 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  the  audtors.  The  present 
report  is  the  third  in  die  series. 
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Structure  of  Ordinary  Ice  Ih 
Part  II:  Defects  in  Ice 
Volume  1:  Point  Defects 

VICTOR  F.  PETRENKO  AND  ROBERT  W.  WHITWORTH 
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INTRODUCTION 

Unlike  most  previous  authors  of  books  on  ice 
physics  (Hobbs  1974,  Fletcher  1970,  Glen  1974),  we 
decided  to  isolate  all  kinds  of  defects  in  ice  into 
separate  reports.  Such  a  division  reflects,  first,  the 
important  role  that  defects  play  in  materials  phys¬ 
ics.  While  most  of  the  physical  properties  of  solids 
are  structurally  sensitive  ones,  in  ice  some  of  them 
cu%  entirely  determined  by  specific  defects.  For 
example,  ion  and  Bjemim  defect  fluxes  are  re¬ 
sponsible  for  electriccil  properties  and  inelastic 
relaxation  in  ice  (see,  e.g.,  Petrenko  and  Ryzhkin 
1984c).  Second,  it  is  in  the  study  of  defects  (dislo¬ 
cations,  point  defects,  grain  boundaries)  that  sig¬ 
nificant  progress  has  been  achieved  since  the  last 
monograph  on  ice  physics  was  published  (Hobbs 
1974).  Because  of  this,  special  attention  will  be 
given  to  the  results  of  studies  from  the  last  two 
decades. 

Following  a  tradition  in  materials  science,  we 
divide  consideration  of  defects  into  three  parts: 
point  defects,  dislocations  and  planar  defects.  In 
this  report  we  discuss  the  structure  and  physical 
properties  of  zero-dimensional  defects,  known  as 
point  defects.  In  the  second  report,  entitled  Struc¬ 
ture  of  Ordinary  Ice  Ih:  Part  II:  Defects  in  Ice;  Volume 
2:  Linear  and  Planar  D^ects,  we  consider  disloca¬ 
tions  that  are  one-dimensional  defects,  and,  fi¬ 
nally,  grain  boimdaries  and  stacking  faults  (two- 
dimensional  defects)  will  be  described.  The  struc¬ 
ture  and  properties  of  another  two-dimensional 
"defect" — the  surface — are  treated  in  a  separate 
report. 

GENERAL  CONSIDERATION 
OF  POINT  DEFECTS 

We  begin  our  review  of  defects  in  ice  with  the 
smallest  and  simplest  ones.  These  are  the  point 
defects  that  are  usually  vacancies,  interstitials, 
impurity  atoms,  bivacancies,  Schottky  and  Frenkel 
pairs,  etc.,  but  in  ice  also  include  special  kinds  of 


defects  associated  with  the  disorder  of  the  hydro¬ 
gen  atoms.  The  term  "point"  does  not,  of  course, 
mean  that  the  defects  do  not  have  a  size,  but  tiKat 
they  are  located  at  a  specific  point  within  the 
crystal.  They  are  usually  locat^  within  one  or 
two  uiut  cells,  although  elastic  deformation  pro¬ 
duced  by  the  defects  can  spread  to  much  larger 
distances. 

Point  defects  in  ice  (as  in  other  solids)  play  a 
fundamental  role  in  many  physical  processes  such 
as  diffusion,  anelastic  relaxation,  electrical  con¬ 
ductivity,  dielectric  permittivity  and  opticed  ab¬ 
sorption.  They  differ  from  one  another  in  struc¬ 
ture,  conc6.itration  and  mobility,  but  in  this  intro¬ 
duction  we  will  consider  some  general  results 
applicable  to  them  all.  These  are  the  relationships 
between  concentration,  formation  energy,  jump 
rates,  diffusion  coefficient  and  mobility. 

Since,  for  a  given  number  of  molecules  the 
arrangement  with  minimum  energy  is  perfect, 
with  no  defects,  some  additional  energy  £f  is  re¬ 
quired  to  produce  a  defect.  For  the  majority  of 
point  defects  in  ice,  £f  is  of  the  order  of  0.5  eV. 
These  values  are  much  greater  than  the  aver¬ 
age  energy  of  thermal  vibration,  ksT  =  0.023  eV 
(-10°C),  so  that  at  first  sight  the  formation  of  these 
defects  is  imjustified.  However,  according  to  ther¬ 
modynamics  at  constant  temperature  and  con¬ 
stant  volume,  thermal  equilibrium  is  attained  at 
minimrim  free  energy  f  rather  than  at  minimum 
internal  energy  U,  where 

F=U-TS.  (1) 

The  appearance  of  defects  increases  the  internal 
energy  U,  but,  because  of  the  large  munber  of 
different  ways  in  which  the  defects  can  be  located 
in  the  crystal,  their  entropy  increases  more  rap¬ 
idly,  and  the  minimum  of  F  occurs  at  a  finite  defect 
concentration. 

We  can  illustrate  this  point  with  a  simple  ex¬ 
ample.  Consider  a  crystal  lattice  consisting  of  N 
identical  points  (atoms,  molecules  or  elementary 
cells).  A  point  defect  with  energy  of  formation  Ef 


could  appear  with  equal  probability  at  any  of  the 
points.  TTien,  the  free  energy  attributable  to  the 
formation  of  n  defects  will  be 

Fd  =  Efn-TSc  (2) 

where  Sc  is  the  configurational  entropy  of  the  n 
defects  in  the  lattice 


Sc  =  itBlnW 


(3) 


and  W  is  the  number  of  arrangements  of  « identical 
defects  over  these  N  identical  sites 


W  = - M! - . 

n!  (N  -  n) ! 


(4) 


At  equilibrium 

^  =  0  (5) 

an 

and  substituting  eq  2-4  into  eq  5  and  using  the 
Stirling  formula — ^ln(N!) = N  InN-N  for  N  » 1 — ^we 
obtain 


— s— =exp 
N-n 


(6) 


Assuming  Ef »  k^T,  we  obtain 


I  kg' 


(7) 


To  be  more  precise,  the  appearance  of  a  point 
defect  not  orUy  changes  the  configurational  en¬ 
tropy,  it  also  ^fects  the  vibrational  entropy  be¬ 
cause  of  a  slight  change  in  the  elastic  lattice  vibra¬ 
tion  spectrum.  This  vibrational  entropy  change  Sf 
would  lead  to  an  additional  factor  of  exp(Se4:B)  in 
eq  7  so  that 


n  =  Wexp 


Calculations  of  Sf  for  defects  are  quite  difficult, 
and  are  far  beyond  this  approach.  We  indicate  here 
only  that,  according  to  experimental  results,  the 
ratio  Sf/k^  may  not  be  small  and  must  be  taken  into 
accovmt. 

The  majority  of  point  defects  in  ice  have  a  high 
mobility  at  "not  too  low"  (T^  100  K)  temperatures, 
migrating  from  one  lattice  site  to  another  by  jump¬ 
ing.  The  equilibrium  positions  of  point  defects  in 
the  lattice  are  separated  by  potential  energy  barri¬ 
ers  of  height  Um,  and  motion  occurs  by  thermally 


activated  hops  over  these  barriers.  Let  us  consider 
the  motion  of  some  electrically  charged  defect 
with  charge  q,  which  we  will  refer  to,  for  brevity,  as 
an  ion. 

The  simple  case  of  the  one-dimensional  motion 
of  such  an  ion  is  shown  in  Figure  1.  To  move  one 
interatomic  distance  in  the  direction  of  d\e  field  E, 
the  positive  ion  would  need  an  activation  energy 
Um-qEa  /2,  and  to  move  against  the  field,  it  would 
need  11^  +  qEa/2.  The  frequency  of  ion  jumps  in 
the  direction  of  the  field  is  expressed  by 


/i  =  V  exp 


kg 


(9) 


and  against  the  field  by 


^2=  vexp 


Um-qaEll 

kg 


(10) 


where  v  is  related  to  the  frequency  of  oscillations 
of  the  ion  in  its  potential  well.  Since  every  jump 
moves  the  ion  over  the  distance  a,  the  mean  drift 
velocity  of  the  ion  in  the  direction  of  the  field 
becomes 

Cd>  -  « (f.  -/2)  -  2v  ,  exp  (-  SW.  j0j . 

(11) 


Since  usually 

qEa«2kBT  (12) 

we  can  substitute  sinh  with  its  argument  in  eq  11 


<Pd>  = 


vqa^ 

kg 


exp  - 


ksT 


E. 


(13) 


Thus,  the  ion  mobility,  which  is  defined  as  the  ratio 
between  drift  velocity  and  electric  field  strength  E, 
will  be 


and  is  independent  of  E.  Um  is  called  the  activation 
energy  of  migration. 

The  quantity 

Vh4v»p(-^j.  (15) 


2 


X 


Figure  1.  Potential  energy  of  an  ion  in  a  crystal  lattice  V/  as  a 
function  of  coordinate  x  vnth  (E>0)  and  without  (E  =  0)  an 
electric  field;  a  is  an  interatomic  distance. 


represents  the  probability  of  jumping  or  the  mean 
nvimber  of  jumps  in  a  given  direction  per  second, 
and  the  dif^ion  coefficient  of  the  ions  is  D= 

For  chaotic  thermal  motion  of  this  kind,  the  mobil¬ 
ity  p  and  the  diffusion  coefficient  D  are  always 
related  by  the  Einstein  relation 


so  that 

D  =  (va2)exp|-|^j.  (17) 

Equation  17,  which  does  not  contain  q,  is  of 
course  also  valid  for  neutral  defects  such  as  vacan¬ 
cies  and  interstitials  and  for  neutral  atoms  and 
molecules  dissolved  in  ice.  When  the  above  simple 
formulae  are  for  some  reason  not  able  to  describe 
the  motion  or  concentration  of  some  specific  de¬ 
fect  in  ice,  we  shall  perform  an  additional  analysis. 

MOLECULAR  DEFECTS 

Tne  principal  molecular  point  defects  in  ice  are 
vacancies  and  interstitials.  A  vacancy  is  a  site  from 
which  a  water  molecule  is  missing,  and  an  intersti¬ 
tial  is  an  extra  molecule  inserted  in  a  space  be¬ 


tween  the  normal  molecules  of  the  structure.  These 
defects  can  clearly  be  formed  in  pairs  by  the  dis¬ 
placement  of  a  molecule  off  its  normal  lattice  site 
to  an  interstitial  position,  but  they  can  also  be 
formed  independently.  For  example,  a  vacancy 
appears  when  we  remove  a  water  molecule  from 
its  regular  position  and  place  it  on  the  surface. 
Dislocations  and  grain  boundaries  can  also  emit  or 
absorb  vacancies  and  interstitials. 

In  thermal  equilibrium  at  a  given  temperature, 
the  defect  concentrations  are  determined  by  their 
formation  energy  Ef  (and  any  associated  entropy 
Sf).  To  a  first  approximation,  the  energy  required 
to  form  a  vacancy  is  equal  to  the  energy  required 
to  remove  a  molecule  from  the  bulk  of  the  crystal 
and  place  it  at  a  kink  on  a  surface  step .  In  so  doing, 
the  number  of  hydrogen  bonds  per  molecule  is 
reduced  from  four  to  two.*  Subsequently,  there 
will  be  an  elastic  relaxation  of  the  lattice  around 
the  vacant  site,  which  slightly  lowers  the  energy. 

The  formation  of  an  interstitial,  on  the  other 
hand,  involves  the  transfer  of  a  molecule  from  a 

•Devlin  (1992a)  recently  found  from  spectroscopic  data  that 
the  surface  molecules  of  crystalline  ice  are  three-coordinate 
bonded,  i.e.,  there  is  only  one,  not  two,  dangling  O-H  bond  per 
molecule.  Nevertheless,  if  a  crystal  is  built  up  by  adding 
molecules  to  the  surface,  the  munber  of  bonds  is  two,  irrespec¬ 
tive  of  the  number  of  dangling  bonds  at  the  surface. 
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(19) 


kink  site  on  the  sutface  to  an  ii  terstitial  position. 
This  again  breaks  two  hydrogen  bonds  and  in¬ 
volves  a  short-range  repulsive  interaction  among 
nearest  neighbors  and  an  elastic  distortion  of  the 
surrounding  lattice.  These  considerations  suggest 
that  theenergy  of  formation  of  these  defects  will  be 
of  about  the  same  magnitude  as  the  cohesive  en¬ 
ergy  per  molecule,  which  is  known  to  be  0.6  eV. 

More  accurate  calculations  using  computer 
simulations  were  made  by  Cotteril  et  al.  (19^).  In 
their  work  they  used  a  realistic  potential  of  inter¬ 
action  between  water  molecules  introduced  by 
Ben-Naim  and  Stillinger  (1972).  Their  calculations 
were  only  restricted  by  their  assumption  that  the 
molecules  were  absolutely  rigid  and  the  small  size 
of  the  molecular  cluster  used  in  their  model.  The 
energies  of  molecular  point  defects  (and  also  of  the 
Bjerrum  L-  and  D-defects  to  be  described  later) 
calculated  by  Cotteril  et  al.  are  svunmarized  in 
Table  1. 

Water  molecules  can  diffuse  through  ice  by 
either  a  vacancy  or  an  interstitial  process.  To  be 
specific,  let  us  consider  such  self-diffusion  as  oc¬ 
curring  by  the  motion  of  vacancies.  The  diffusion 
flux  of  vacancies /v  is  given  by  the  formula 

;v  =  -Dvgrad(«v)  (18) 

where  Dy  is  the  diffusion  coefficient  of  vacand^ 
defined  in  eq  17  and  tty  is  the  vacancy  concentra¬ 
tion.  A  particular  water  molecule  in  the  ice  can 
migrate  by  vacancy  diffusion  only  if  one  of  its 
nearest-neighbor  lattice  sites  happens  to  be  va¬ 
cant.  Consequently,  the  coeffident  of  self-diffu¬ 
sion  of  molecules  Ds  is  much  smaller  than  the 
coeffident  of  vacancy  diffusion 


Table  1.  Defect  formation  energy  (e  V)  in  ice  (after 
Cotteril  etal.  1973). 


Cluster  size  and  elastic  relaxation 


D^t 

16,  unrelaxed 

16,  relaxed 

432,  unrelaxed 

Cohesive 

0.466 

— 

0.465 

energy 

Vacancy 

0.466 

0.447 

0.465 

L-defect 

0.467 

0.388 

0.403 

D-defect 

0.491 

0.407 

0.428 

Interstitial 

0.610 

0.528 

0.695 

In  practice,  the  self-diffusion  coeffident  of  a  mate¬ 
rial  can  be  determined  using  a  tracer  technique, 
provided  that  the  material  possesses  suitable,  usu¬ 
ally  radioactive,  isotopes.  For  ice,  three  isotopes 
have  been  used  in  sudi  difiusion  measurements: 
deuterium,  ^H;  tritium,  ^H;  and  i®0.  The  general 
procedure  is  to  place  a  known  concentration  of  the 
isotope  onto  one  end  of  a  block  of  ice,  which  is  then 
left  at  a  fixed  temperature  for  a  certain  period  of 
time.  The  distribution  of  the  isotope  along  a  par¬ 
ticular  axis  is  then  measured  experimentally.  Table 
2  and  Figure  2  summarize  the  results  from  such 
experiments. 

If  we  substitute  eq  8  and  17  into  eq  19,  we  obtain 

D..D.exp(-j&)  (20, 

in  which 

D,o=(vfl2)exp(^  (21) 

and 

Es  =  Ef+lJm.  (22) 

Several  trends  in  the  experimental  data  in  Table  2 
are  of  interest.  First,  the  difiusion  coefiidents  of 
oxygen  and  hydrogen  isotopes,  as  well  as  their 
adivation  energies,  are  approximately  equal.  This 
is  major  evidence  that  in  ice  oxygen  and  hydrogen 
diffu^  cooperatively  by  the  diffusion  of  whole 
water  molecules.  The  scatter  in  Ds  and  Es  (rela¬ 
tively  smaU)  could  be  caused  by  different  degrees 
of  imperfection  of  the  ice  used.  For  the  majority  of 
solids,  it  has  been  established  that  the  diffusion 
rate  along  grain  boundaries  or  along  dislocation 
lines  considerably  exceeds  the  bulk  rate.  Hence, 
one  could  expect  larger  values  of  Ds  in  polycrystal¬ 
line  ice  and  in  ice  with  high  dislocation  density 
(which  was  unknown  in  the  papers  cited). 

The  key  question  in  the  study  of  self-diffusion 
in  ice  is  about  whether  the  process  occurs  by  the 
motion  of  vacancies  or  interstitials.  Since,  in  the 
majority  of  crystals  vacandesdominate  interstitials, 
self-diffusion  in  ice  was  for  a  long  time  also  ex¬ 
plained  in  terms  of  vacancy  diffusion.  However, 
experiments  by  Higashi,  Hondoh  and  others  dur¬ 
ing  last  decade  have  clearly  proven  that  the  domi¬ 
nant  process  in  diffusion  in  ice  is  that  ascribable  to 
interstitials.  Their  experiments  involved  the  care¬ 
ful  study  of  the  growth  and  annealing  of  disloca¬ 
tion  loops  produced  by  rapid  cooling. 
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Table  2.  Self-diffusion  coefficient  O*  and  its  activation  energy  £,  in  ice. 


Measurement 

Self  diffusion 
coefficient  Dg 
(m^s-t)at 
-10°C 

Activation 
energy  Eg  (eV> 

Reference 

Oxygen-18  in  polycrystalline  ice. 

1.0  X  10-w 
(at-1.8°C) 

— 

Kuhn  and  Thiirkauf  (1958) 

Deuterium  in  polycrystalline  ice. 

1.0  X  10-w 
(at-1.8“C) 

— 

Kuhn  and  Thiirkauf  (1958) 

Tritium  perpendicular  to  c-axis  in  artificial  single 
crystals. 

1.53  X  Ifi-ts 

0.58  ±0.05 

Dengel  and  Riehl  (1963) 

Tritium  perpendicular  and  parallel  to  c-axis  in 
natural  single  crystals.  Diffusion  time  of 

7  days. 

2.45x10-15 

0.67  ±0.09 

Itagaki  (1964) 

Tritium  parallel  to  c-axis  in  rutural  single 
crystals.  Diffusion  time  greater  than  2  days. 

2.12x10-15 

0.54  ±  0.09 

Itagaki  (1966, 1%7) 

Tritium  perpendicular  to  c-axis  in  natural  single 
crystals.  Diffusion  time  greater  than  2  days. 

2.42  X 10-15 

0.65  ±0.05 

Itagaki  (1966, 1%7) 

Tritium  parallel  to  c-axis  in  artificial  single  crystals. 

2.15x10-15 

0.63  ±0.03 

BlicksetaL(1966) 

Tritium  perpendicular  to  c-axis  in  artificial  single 
crystals. 

2.40x10-15 

0.63  ±0.03 

BlicksetaL(1966) 

Oxygen-18  in  artificial  single  crystals. 

3.1  X 10-15 

0.68  ±0.13 

Delibaltas  et  al.  (1966) 

Tritium  parallel  to  c-axis  in  artificial  single  crystals. 

2.5x10-15 

0.63  ±0.05 

Dengel  et  al.  (1966) 

Tritium  parallel  to  c-axis  in  artificial  single  crystals. 

1.44x10-15 

0.62  ±0.03 

Ramseier  (1967) 

Tritium  perpendicular  to  c-axis  in  artificial  single 
crystals. 

1.68x10-15 

0.62  ±0.07 

Ramseier  (1967) 

Tritium  parallel  to  c-axis  in  natural  single  crystals. 

1.54  X 10-15 

0.61  ±  0.02 

Ramseier  (1%7) 

Tritium  perpendicular  to  c-axis  in  natural 
single  crystals. 

1.65  X 10-15 

0.62  ±  0.04 

Ramseier  (1%7) 

Consider  what  happens  in  the  crystal  during  Dislocation  loops  (Fig.  3)  can  expand  or  con- 

rapid  cooling.  The  term  "rapid"  means  that  the  tract  because  of  diffusion  of  point  defects  toward 

cooling  takes  place  during  a  period  much  less  than  (or  from)  the  loops,  caused  by  a  change  in  temper- 

that  required  for  the  new  equilibrium  defect  con-  ature  or  applied  pressure.  Monitoring  the  growth 

centrations  to  be  established.  When  the  tempera-  or  contraction  of  dislocation  loops  by  means  of  X- 

ture  is  lowered,  equal  numbers  of  vacancies  and  ray  topography  in  ice  monocrystals,  Higashi  and 

interstitials  disappear  rapidly  owing  to  their  mu-  Hondoh  with  co-authors  detenrxined  the  diffusion 

tual  annihilation.  Then  some  fraction  of  defects  coefficient  of  interstitials  in  ice  Dj,  the  formation 

that  are  predominant  must  "die"  in  a  slower  and  energy  of  interstitials  Has,  the  activation  energy  of 

more  difficult  way— for  instance,  by  diffusing  to-  their  motion  Um,  the  activation  volumes  of  inter¬ 
wards  the  surface  or  concentrating  into  disk  clus-  stitials  Vy  and  of  their  migration  Vim,  and  absolute 

ters,  as  shown  in  Figure  3.  The  edges  of  these  disks  values  of  the  interstitial  concentration  nj,  and  the 

are  edge  dislocation  loops,  with  the  Burgers  vector  entropy  Sy/ks.  Table  3  contains  all  of  these  data, 

normal  to  the  plane  of  the  disk.  If  loops  of  a-type  According  to  these  results,  the  activation  en- 

are  predominant,  then  the  major  defect^  are  vacan-  ergy  of  self-diffusion  in  ice  Eg  =  Eas  +  Eim  =  0.56  eV, 

cies;  if  loops  of  b-type  are  predominant,  then  the  and  the  absolute  value  of  the  self-diffusion  coeffi- 

major  defects  are  interstitials.  cient  is  in  good  agreement  with  results  obtained 
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TEMPERATURE  (“C) 


Figure  2.  Comparison  of  the  self-diffusion  co¬ 
efficients  D  in  ice  obtained  by  the  tracer  tech¬ 
nique  (^H  and  and  calculated  from  the 

di^sion  coefficients  of  self-interstitials  ob¬ 
tained  by  the  X-ray  topography  (solid  line) 
(after  Oguro  et  al.  1988). 


IC^fT  (K-‘) 


Figure  3.  Two  types  of  the  prismatic  dislocation  loop — 
(a)  xxicancy  type  and  (b)  interstitial  type.  Schematic 
atomic  arrangement  at  the  top  corresponds  to  edge 
dislocation  loops  of  different  signs  (af^  Oguro  et  al. 
1988). 
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Table  3.  Sdf-diffusion  intentitial  parameten  in  ice  at  (TC 
(after  Hondoh  1992).  Vmol »  a  molecular  volume  in  ice 
(V«oi- 3.3x10-23  cm3). 

tu/N  Em(eV)  SJi(M  Di(nfi/s)  U.feW  V^/V«<  Vim/Wmci 
2.8xl(H  0.40  4.9  2x10^  0.16  -0.26  0.40 


earlier  by  the  radioactive  tracer  technique  (Fig.  2). 
We  refer  readers  interested  in  the  details  of  an 
elegant  experimental  technique,  employed  by 
Higashi  and  Hondoh  et  aL,  to  their  review  papers 
(Hondoh  et  al.  1987,  Oguro  et  al.  1988,  Hondoh 
1992),  where  one  can  find  references  to  their  origi¬ 
nal  papers. 

Thus,  from  measurements  of  self-diffusion  in 
ice  it  is  not  possible  to  determine  any  of  the  char¬ 
acteristics  of  vacancies,  since  the  predominant 
defects  are  interstitials  (at  least  at  T  S  -50°C).  The 
only  studies  in  which  activation  parameters  of 
vacancies  in  ice  were  possibly  determined  were 
those  in  which  a  positron  anrfihilation  technique 
was  used  (MogensenandEldrup  1978,  Eldrup  and 
Morgensenl978).Wheniceisexposedtopositrons, 
thermalized  positrons  and  positronium  atoms 
(electron  +  positron)  are  captured  mainly  by  va¬ 
cancies.  Therefore,  Aeir  lifetimes  depend  on  the 
concentration  and  mobility  of  vacancies.  It  is  ex¬ 
tremely  difficult  to  estimate  the  absolute  accuracy 
of  this  method  since  it  is  based  on  a  considerable 
number  of  assumptions  and  approximations.  The 
vacancy  migration  energy  foimd  in  these  papers 
was  0.^  ±  0.07  eV  and  its  formation  energy  was 
foimd  to  be  between  0.2  and  0.35  eV. 


PROTONIC  DEFECTS 

As  explained  in  the  earlier  reports  on  the  Struc- 
tureof Ordinary  IcelhiPart  I:  Ideal  Stnicture  o//ce  and 
Electrical  Properties  of  Ice  (Petrenko  1993a,b),  many 
importantproperties  arise  fi'om  the  disorder  of  the 
protons  in  the  structure.  This  allows  molecules  to 
rotate,  giving  rise,  for  example,  to  dielectric  polar¬ 
ization,  but  such  rotations  cannot  occur  with¬ 
out  the  aid  of  special  kinds  of  point  defects  in 
die  protonic  suteystem.  The  structure  of  ice  is 
shown  in  Figure  4a.  There  are  two  possible  posi¬ 
tions  for  protons  on  each  bond,  but  in  the  perfect 
structure  the  occupation  of  these  sites  is  governed 
by  the  so-called  Bemal-Fowler  rules  or  simply  the 
"ice  rules": 


1.  Two  protons  are  located  near  each  oxygen 
atom. 

2.  There  is  one  proton  per  hydrogen  bond. 

The  first  rule  ensures  electroneutrality  of  the  water 
molecules,  and  the  second  ensures  the  perf^on 
of  the  hydrogen  bonds.  Violations  of  these  rules  at 
particular  bonds  or  molecules  result  in  the  forma¬ 
tion  of  the  four  protonic  defects  shown  in  Figtue 
4b.  These  are  the  H3O  and  OH~  ions  and  the 
Bjerrum  L-  and  D-defects,  which  are  bonds  with 
no  protons  or  two  protons  respectively  (Bjerrum 
1951).  An  ion  moves  from  one  site  to  another 
merely  by  the  transfer  of  a  proton  from  one  end  of 
a  bond  to  another,  and  a  Bjerrum  defect  moves  by 
the  rotation  of  a  water  molecule. 

The  structure  of  die  hydrogen  bonds,  along 
which  the  protons  move,  constitutes  a  three-di¬ 
mensional  grid.  The  transfer  along  this  grid  of  any 
of  the  above  defects  reorients  water  molecules  in 
such  a  fashion  that  a  second  defect  of  the  same  type 
cannot  pass  along  the  same  path  in  the  same 
direction. 

To  better  understand  what  happens  in  the  pro¬ 
cess  of  ionic  and  ^errum  defect  transfer,  let  us 
consider  the  track  of  an  arbitrary  charge  carrier,  a 
positive  H3O+  ion,  for  instance.  When  we  are 
speaking  about  ionic  or  Bjerrum  defect  transfer,  it 
is  important  to  realize  that  oxygen  atoms  do  not 
"hop"  from  one  lattice  site  to  the  next  one.  Only 
protons  move,  widi  a  subsequent  redistribution  of 
the  electronic  density  on  the  neighboring  mol¬ 
ecules.  Thus,  when  speaking  about  HsQ*'  motion, 
we  have  in  mind  that  a  characteristic  distribution 
of  three  protons  near  an  oxygen  atom  shifts  in 
space. 

When  a  proton  jumps  along  a  hydrogen  bond 
from  one  water  molecule  to  anoAer,  a  pair  of 
HaO  and  OH-  ions  is  produced,  as  shown  in 
Figure  4b.  Formation  of  this  ion  pair  and  further 
separation  of  the  ions  require  considerable  energy 
of  about  1  eV  (see  the  review  of  experimental 
results  in  Petrenko  1993a).  After  an  ion  is  formed, 
it  can  migrate  from  one  site  to  another  through 
proton  jumping  along  hydrogen  bonds.  This  pro- 


b.  Ice  rvith  HiO*  and  OH-  ions  and  D-  and  L-defects. 
Figure  4.  Ice  structures. 


cess  is  illustrated  schematically  for  the  straight¬ 
ened  out  chain  of  molecules  shown  in  Figures  5a 
andb. 

When  a  proton  jiunps  from  one  bond  to  <m- 
other,  two  neighboring  defective  hydrogen  bonds 
are  produced.  These  are  the  so-called  D-  and  L- 
defects.  Since  a  proton  has  a  positive  electrical 
charge,  it  is  not  difficult  to  see  that  a  D-defect  is 
positively  charged  and  an  L-defect  is  negatively 
charged.  Formation  and  separation  of  a  pair  of 
Bjerrum  defects  require  energy  of  0.68  eV.  Already 
formed  L-  and  D-defects  can  migrate  in  the  ice 
lattice  by  protons  jumping  between  neighboring 


hydrogen  bonds  as  illustrated  for  a  D-defect  in 
Figures  5c  and  d. 

An  H3O+  ion  migrates  from  left  to  right  along  a 
water  molecule  chain  by  consecutive  jumps  of 
protons  along  hydrogen  ^nds.  Therefore,  if  in  the 
beginning  the  protons  on  all  bonds  in  the  chain 
were  located  in  the  left-hand  positions  (Fig.  5a), 
after  an  ion  has  passed,  all  protons  will  be  in 

the  right-hand  positions.  This  path  is  now  closed 
for  proton  transfer  in  the  same  direction.  The 
reader  can  easily  make  sure  that  the  migration  of 
an  OH-  ion  from  right  to  the  left  would  result  in 
exactly  the  same  rearrangement  of  this  molecular 


b.  Cham  after  the  ion  has  passed.  Obviously,  no  fiirther  H3O  ion  can  go  along  the  same 
chain  in  the  same  direction. 


c.  D-defect  motion  from  to  right  along  the  chain  as  in  Figure  5b.  The  arrows  indicate 

successive  hops  of  protons  from  one  bond  to  next  one. 


d.  Chain  after  the  D-d^ect  has  passed.  Comparison  with  Figure  5a  shows  that  motion  pfD- 
d^ect  restored  positions  of  the  protons  along  the  chain.  It  is  ready  (or  "open")  now  for  the 
passage  of  another  ion. 

Figure  5.  Schematic  representation  ofHsO*-  ion  (a,  b)  and  D-d^ect  (c,  d)  motion  along  a 
chain  of  hydrogen  bonds  in  ice. 


chain.  Since  H3O+  and  OH-  ions  drift  in  opposite 
directions  in  an  external  electric  field,  we  come  to 
the  conclusion  that  both  types  of  ions  drifting  in  a 
given  field  "block"  molecular  tracks  in  the  same 
maimer. 

If  we  now  consider  the  migration  in  the  same 
field  of  a  positively  charged  EMefect  along  the 
same  chain  (Fig.  5c)  "blocked"  by  an  HaO*  ion 
passing  along  it,  it  turns  out  that  the  D-defect 
"unblocks"  the  chain,  restoring  the  initial  configu¬ 
ration  to  that  of  Figure  5d,  which  is  the  same  as 
Figure  5a.  In  the  same  fashion  the  migration  of  a 


negatively  charged  L-defect  "imblocks"  the  chain 
if  it  migrates  from  right  to  left. 

Ions  and  ^errum  defects  migrating  through  ice 
reorient  water  molecules,  altering  conditions  for 
their  transfer.  One  can  say  that  in  ice  protonic 
defects  interact  in  a  very  special  way:  via  those 
changes  in  the  molecule  configuration  that  they 
produce  during  their  transfer.  In  addition  the  more 
usual  interactions  between  defects  also  take  place. 
First,  protonic  defects  are  electrically  chuged, 
which  results  in  a  strong  Coulomb  interaction,  and 
second,  as  we  shall  see  later  (in  the  Activation 


Volume  of  Protonic  Defects  section),  die  presence  of 
a  defect  introduces  a  lattice  defonnaticHi  that  can 
lead  to  interactions  via  elastic  stresses. 

In  thermal  equilibrium,  which  is  characterized 
by  maximum  oitropy  S,  water  molecules  are  ori¬ 
ented  randomly,  subject  to  the  ice  rules.  Any  flux 
of  defecte  of  the  same  type  (H3O  ions,  for  ex¬ 
ample)  ;  I  will  increase  Ae  number  of  molecules 
oriented  along  the  direction  of  the  flax.  The  new 
configuration  of  water  molecules  is  more  ordered 
and  ^erefore  has  smaller  entropy.  These  ideas 
were  formulated  mathematically  by  Jaccard  (1959, 
1964)  and  independently  by  Otisager  and  Depuis 
(1960).  They  introduced  the  configurational  vector 
Q, describing  changes  occurring  In  ice  because  of 
protonic  defect  transfer.  If  at  time  f  =  0  ice  is 
impolarized  and  in  thermal  equilibrium,  then 


^  =  (/I  -h-h  +/4 

Jo 


where  (i  =  1,2^4)  denote  the  fluxes  of  H3O+, 
OH-  and  D-  and  L-defects  respectively.  This  nota¬ 
tion  of  protonic  defects  is  generally  accepted  and 
we  shall  use  it  hereafter.  If  we  introiduce  die  quan¬ 
tities  iii  given  by 

tl,  =  1,-1, -1,1  for  i  =  l,2A4,  (24) 

then  eq  23  can  be  simplified  to 


-*  S  r‘ 

i=iio 


The  choice  of  +  or  -  before  in  eq  23  accounts 
for  different  effects  produced  by  defects  on  die 
reorientation  of  water  molecules.  Thus,  for  ex¬ 
ample,  as  we  have  seen  in  Figure  5,  the  effect  of  an 
HaO^  ion  (i  =  1)  is  opposite  to  the  effect  of  a  D- 
defect  (i  =  3). 

Configurational  entropy  Sc  must  somehow  be 
related  toll  .For  small  decrements  of  entropy  we 
can  use  the  first  term  of  a  Taylor  series 


The  proporticmality  foctor  in  eq  27,  which  we 
doiote  as  ^/T,  can  in  principle  be  both  calculated 
dieoretically  and  detamined  experimentally.  At¬ 
tempts  to  cakulate  ^  theoretically  Jaccard  1959, 
1964;  Hubmann  1979a,b)  yieldonly  die  ri^t  order 
of  magnitude  oi  diis  quantity,  and  die  most  reli¬ 
able  value  is  that  found  in  an  analysis  of  experi¬ 
mental  data  by  Hubmatm  (1979a) 

|  =  3.85kBrco  (28) 

where  Too  =  2.76  A  is  the  oxygen-oxygen  distance 
in  ice. 

According  to  Ryzhldn*  the  difierence  between 
the  theoretical  and  experimental  values  of  O/T  is 
caused  by  taking  into  account  reorientaticHi  of 
molecules,  not  <xily  in  this  water  molecule  chain, 
but  also  reorientation  of  nei^iboring  molecules 
when  considering  a  defect  transfer  along  a  chain. 
However,  sudi  calculations  are  quite  cumbersome 
and  have  not  been  made  yet.  Accordingly,  we  shall 
use  the  value  of  ^  found  experimentally. 

Charge  carriers  in  ordinary  ice  ht  oie,  under 
usual  conditions,  one  of  these  four  types  of  pro¬ 
tonic  defects.  A  detailed  discussicai  on  die  nahue 
of  electric  charges  in  ice  can  be  fotmd  in  die  report 
by  Petroiko  (1993a).  "Usual  ccmditions"  are  un¬ 
derstood  here  as  an  absence  of  excitation  by  direct 
ultraviolet  light,  electron  beanig,  y  rays,  etc. 

The  electric  current  density/  inicecannowbe 
written 

T  =  i  e;  A  <29) 

i  =  l 

where  e;  are  the  efiective  charges  carried  by  the 
defects.  Since  HjO^  and  OH-  ions  and  L-  and  D- 
defects  are  created  and  annihilated  in  pairs,  it  is 
obvious  that 

ei  =  -e2  and  C3  =  -e4 .  (30) 

If  we  place  an  extra  "foreign"  proton  into  an  odi- 
erwise  normal  ice  lattice,  we  will  inevitably  pro¬ 
duce  two  carriers  at  once:  an  H3O*  ion  and  a  D- 
defect,  as  shown  in  Figure  6.  Hence,  the  sum  of 
these  charges  equals  a  proton  charge  e 


There  is  no  linear  term  because  ASc  cannot 
depend  on  the  sign  of  the  fluxes  producing  it.  For 
the  time  derivative  of  the  entropy  we  have 

S«-nQ.  (27) 


ei  +e3  =  e. 


*Personal  communication  widt  Dr.  Ivan  Ryzhkin  during  dw 
International  Symposium  on  Physics  and  Chemistry  of  ke^ 
September  1991,  Sapporo,  Japan. 
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Figures.  Crettti(m<iftxvopomtdefeds(H30*’ionitndD-defect)whertaprotonisplttced 
into  defect-free  ice  lattice. 


To  find  out  how  the  elementary  protonic  charge  is 
shared  between  the  ion  and  the  D-de£ect,  we  must 
consider  both  the  fraction  of  the  intermolecular 
distance  that  the  proton  travels  in  the  case  of  ionic 
or  ^errum  defect  transfer  and  die  associated  die- 
placement  of  electronic  charge  density  when  the 
defect  moves  along  the  chain. 

The  most  strict  theoretical  analysis  of  this  prob¬ 
lem  was  conducted  by  Scheiner  and  Nagle  (19^) 
who  foimd 

e3  =  (036±0.03)e.  (32) 

The  other  way  to  find  ei  andesis  by  experiment. 
In  ice  the  charge  of  the  predominant  defect  is 
incorporated  in  the  formula  for  static  dielectric 
permittivity  Es  (see  Petrenko 1993a).  Using  this  fact 
and  analysis  of  experimental  data  of  measure¬ 
ments  of  £s  for  specially  doped  ices  (in  such  a 
fashion  that  either  ions  or  ^errum  defects  pre¬ 
dominated),  Hubmaim  (1979a,b)  found 

e3  =  (0.38±0.01)c  (33) 

which  is  remarkably  close  to  the  calculations  of 

Scheiner  and  Nagle.  Other  recent  experimental 

results  for  ei  and  ea  are  given  in  Table  4.  These 
results  are  close  to  those  represented  in  eq  32 
and  33. 


The  rate  of  heat  generation  in  a  unit  volume 
Q  by  current  flowing  throu^  ice  is  determined 
by,  two  efiects.  The  first  is  simply  ttie  Joule  heat 
£/,  and  the  second  is  the  change  in  the  configura- 


tional  entropy  TSq  Accordingly 

Q  =  Ef-<tfaa. 

(34) 

Using  eq  25  and  29  we  obtain 

Q=  i  (e,E-TijC»n)7,. 

(35) 

i=l 


As  canbe  seen  from  eq35,  the  efiiective  force  acting 
upon  cltarge  carriers  in  ice  consists  of  two  terms,  as 
opposed  to  ordinary  conductors 


Table  4.  Ratios  of  effective  electric  charge  of 
HaO  ion  (ei)  and  of  D-def^  (ea)  to  proton 


charge  (e). 

Ci/e 

Reference 

0.61  ±0.02 

0.39  ±0.01 

Zarelskii  et  aL  (1988) 

0.62  ±0.01 

0.38  ±0.01 

Hubmann  (1979a) 

0.626 

0.374 

Take!  and  Maeno  (1987) 

0.64  ±0.03 

0.36  ±0.03 

Scheiner  and  Nagle  (19K) 

(41) 


Fi=eiE-n,<l>Q.  (36) 

Accordingly,  for  the  charge  carrier  fluxes  we  have 

j  i  =  (e,  E-  ti,4>  Q  )  n,  /«,  (37) 

where  (li  is  the  mobility  (i.e.,  velocity  per  uiut 
electric  field)  of  defect  i.  If  the  distribution  of 
protonic  defats  in  ice  is  not  uniform,  we  have  to 
add  the  diffusion  flux  to  eq  37,  giving 

y,  =  (^,£-Ti,4>n)  a,  p,/e,-D,grad  n,  (38) 

in  which  the  diffusion  coefficient  Di  is  equal  to 
itsTV/ki  I-  '11^6  lour  equations  (eq  38,  i  =  1, 2, 3, 4) 
are  the  most  general  equations  describing  charge 
carrier  transfer  in  ice  and  can  account  for  all  elec¬ 
trical  properties  of  ice  known  today  (Petrenko 
1993a). 

Statistics  of  protonic  point  defects  in  ice 

To  calculate  the  equilibrium  ccmcentration  of 
protonic  point  defects,  we  can  follow  the  same 
general  procedure  that  was  described  in  the  Gen¬ 
eral  Consideration  of  Point  Defects  section.  For  ice 
with  ionic  defects,  such  a  calculation  was  made  by 
Ryzhkin  (1985),  who  carefully  calculated  the  en¬ 
tropy  of  ice  with  H3O  and  OH~  ions.  Minimizing 
the  free  energy  F,  Ryzhkin  then  obtained  for  the 
equilibrium  concentration  of  ions 

=2.exp(— (39) 

N-Tni  3  \  2*^/ 

where  £»  is  the  formation  energy  of  a  pair  of  ions, 
and  since  Eai » 

As  we  can  see,  the  specific  structure  of  ice  mani¬ 
fested  itself  only  in  the  appearance  in  eq  40  of  a 
factor  of  2/3  before  the  exponent.  We  can  interpret 
this  fiictor  in  simple  terms  as  the  ratio  of  the 
number  of  possible  orientations  of  an  ion  (four 
available)  to  the  number  of  possible  orientations  of 
a  water  molecule  (which  is  six).  The  factor  1  /2  is  in 
the  exponent  because  ions  are  created  in  pairs. 

Ihe  general  case  of  the  presence  of  ions  and 
^errum  defects  was  also  considered  by  Ryzhkin 
(impublished  data,  1985).  Let  us  denote  the  !> 
def^  concentration  as  no  and  the  concentration 
of  Lrdefects  as  nt.  If  D-  and  L-defects  ate  formed 
by  thermal  fluctuations,  i.e.,  in  pairs,  then  obvi¬ 
ously 


«D  =  «L- 

Using  a  method  similar  to  the  one  seen  above, 
Ryzhkin  found  that  in  the  presence  of  both  types  of 
defects 

S/(kBN)  =  -2x  ln(x)  -  (1-2  x)  ln(2  (1-2  x)/3) 
-4yln(2y)-2(l-2y)ln(l-2y)  (42) 

where 

X  =  nfN  and  y  =  no/N .  (43) 

Minimizing  the  free  energy  of  ice  with  the  entropy 
from  eq  42,  we  obtain  for  the  ^errum  defect  con¬ 
centration 

nD=  "L*  "H20  e*P  (“^^)  • 

Protonic  defects  can  also  be  introduced  into  ice 
by  doping  with  certain  impurities.  This  case  will 
be  considered  in  the  Impurities  section. 

Electrical  properties  of  protonic  defects 

As  we  have  indicated  previously,  protonic  de¬ 
fects  are  charge  carriers  in  ice  and  hence  determine 
its  electrical  properties.  A  detailed  description  of 
file  electrical  properties  of  ice  and  experimental 
techniques  for  their  investigation  can  be  found  in 
the  recent  review  by  Petrenko  (1993a).  Here,  we 
shaU  only  point  out  what  information  concerning 
proUmic  defects  can  be  obtained  from  aiudysis  of 
the  electrical  properties  of  ice. 

When  an  alternating  electric  field  with  angular 
frequency  a>  is  applied  to  ice,  in  addition  to  the 
drift  current  (eq  29^  J^,  we  have  to  add  the  dis¬ 
placement  current  so  that  the  total  current  is 

/  =  y  dis  +  *  (45) 

/dis  =  C0£«-^  (46) 

where  co  is  the  dielectric  permittivity  of  free  space 
and  e.  is  the  high-frequency  dielectric  permittiv¬ 
ity  of  ice  arising  firom  its  electronic  polarizability; 
&.s3.2. 

A  complete  description  of  the  electrical  proper¬ 
ties  of  ice,  i.e.,  its  electrical  conductivity  a((o)  and 
dielectric  permittivity  e((o),  can  be  foimd  in  the 
system  of  eq 25, 29, 45  and  46  (see  Petrenko  1993a, 
for  instance).  In  the  absence  of  interfaces  and  con¬ 
centration  gradients,  the  answer  can  be  expressed 
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as  follows 


(55) 


<T(a))  =  <^  + 


(<y- -<%)(«) 

1  +  (<o  xd)2 


e  (w)  =  + 


«s-«- 
1  +  (<0  Xd)^ 


(47) 

(48) 


Here,  we  have  introduced  new  terms 

4 

=  X  Oi  with  oj  =  rti  |ej  n, 


i  =  l 


e2/flr.  = 


Oi  +  05  <^  +  04 


=  0)0  =  ® 


(49) 

(50) 

(51) 


£s  =  £-  + 


oi  +  <^  <^  +  04 


^3 


£o4» 


0l  +  +  04)^ 

ei2  ^32  ) 


(52) 


By  varying  the  temperature  or  the  type  and 
concentration  of  dopants,  or  botii  (see  the  Impuri¬ 
ties  section),  any  one  of  the  four  terms  can  be 
made  dominant.  Then  eq  52  is  simplified  to  yield 

el 

es=£-  +  -^— •  (53) 

Co® 

Measurements  of  £$  for  doped  ice  in  which  iotuc 
defects  and  Bjerrum  defects  dominate  the  electri¬ 
cal  properties  thus  permit  both  the  charges  Ci  and 
^  to  be  determined.  This  is  the  basis  of  the  results 
quoted  in  Table  4. 

As  we  can  see  firom  eq  49  and  50,  ou  is  deter¬ 
mined  by  the  largest  Oi,  and  0$  is  determined  by 
the  lower  of  (oi  +  and  (o’a  +  04)  .  In  pure 
monocrystalline  ice  (c^ + 04) »  (oi  +  05)  and  a\  »  c^. 
Consequently,  the  measurement  of  a.,  as  a  func¬ 
tion  of  temperature  allows  us  to  determine  the 
activation  energy  of  the  partial  conductivity  04  of 
L-defects,  E04 

O'-  “  04  =  <^40  exp  (-  .  (54) 

Since  the  conductivity  of  Bjerrum  defects  greatly 
exceeds  (usually  by  two  or  three  orders  of  magni¬ 
tude)  the  ionic  conductivity,  and  the  mobility  of 
HaO^  ions  exceeds  the  mobility  of  OH~  ions,  we 
have 


■'-(irf"' 

Iherefore,  from  fite  temperature  dependence  o/D, 
it  is  possible  to  calculate  the  ccmductance  activa¬ 
tion  energy  of  HaO*  ums.  Since  omductivities 
are  the  pr^ucts  of  ccvioentrations  and  mobilities, 
the  activatitm  energies  £0,  obtained  are,  generally 
speaking,  the  siuiiis  of  formation  energies  (Eti/2) 
and  activation  energies  of  migration  of  ttiese  de¬ 
fects  Enu.  To  separate  out  ( Eaj/2)aiKl  Em,  usually  a 
known  concentration  of^)edfic  defects  is  achieved 
by  some  technique  and  fiten  die  are  measured. 
The  known  concentration  of  dte  defects  can  be 
created  wid\proton  injection  (PetrenkoetaL  1983), 
afield  efiecttransistor(Petra\koandMaeno  1987) 
or  with  doping  of  ice  (see  tl^  Impurities  sectkm). 
Thus,  in  many  cases  it  is  possible  to  determine 
separately  concentratimis  n,-  and  mobilities  p,-  of 
protoxuc  defects. 

The  most  well-established  omductivity  param¬ 
eters  of  pure  monocrystalline  ice  are  ^  high- 
frequency  conductivity  ou  and  its  activation  en¬ 
ergy  Eo-.  For  T  >  ^*C  these  depend  oiily  mi  tihe 
L-defects.  These  quantities  canbe  well  reproduced 
because,  first,  high-fiequency  measurements  do 
not  dep>end  upMi  fire  t)^  of  electrode  used  and, 
second,becauseofthecomparativeeasewithvshkh 
specimens  can  be  obtain^  that  are  pure  enou^ 
for  the  L-defect  concentration  to  be  intrinsic.  The 
same  holds  for  Xd  and  its  activation  energy  Et, 
sirKe  according  to  eq  49-^51,  Xd~2  oe  ou  when  o-  » 
(jg  ( i.e.,  when  the  majority  type  of  charge  carrier  is 
ctearly  distinguished).  These  four  values  at  T  = 
-lO'C  are  as  follows 

Eo-=  (0.58±0.03)eV 

Er  =  (058±0.01)eV 

Oi.  =  (1.8  ±0.2)10-5  £2-1  m-l 

XD  *  5xl0-5s.  (56) 

The  Debye  dispersion  shown  in  Figure  7  is  repro¬ 
duced  very  well  in  good  quality  single  crystals: 
pure,  homogeneous  and  annealed  (for  relaxation 
of  internal  stresses). 

Much  greater  scatter  is  observed  in  measure¬ 
ments  of  the  static  conductivity  '  mdoped  ice 

and  its  activation  energy  Eo^.  Values  of  Ea^  have 
been  foimd  to  lie  between  0  and  0.7  eV  and 
(-lO'C)  varies  firom  10-^  down  to  6  x  lO-i®  £2"i 
m-i .  This  variation  is  explained  by  the  difficulty  of 
purifying  ice  to  a  level  at  which  intrinsic  ions 
dominate  over  ions  supplied  by  impurities.  That  is 
why  von  Hippel  et  al.  (1973)  concluded  that  there 
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Figure  7.  Typical  frequency  dependence  of  permittivity  t  and  cf  conductivity  a  of 
pure  single  crystalline  ice  atT  =  -10C°(t^ter  Petrenko  1993a). 


are  no  intrinsic  ions  in  ice  and  that  all  ions  origi¬ 
nate  from  impurities.  If  we  refer  to  Table  I  of  their 
paper  (von  Hippel  et  al.  1971)  for  the  value  of 
(2.5  X  1(H  m"i)  and  compare  it  with  the 
minimum  observed  value  of  6.4  x  10-1°  mr\ 
then  we  might  agree  with  diem  that  their  ions  were 
of  impurity  origin.  The  question,  however,  re¬ 
mains  open:  were  the  specimens  in  the  experi¬ 
ments  of  Worz  and  Cole  (1969)  or  Petrenko  et  al. 
(1983)  pure  enough  that  their  conductivity  may  be 
considered  intrinsic?  We  observed  the  following 
general  trend:  the  purer  the  ice  is,  the  smaller  is 
and  the  greater  is  Eas.  Therefore,  the  value  of  the 
activation  energy  of  creation  of  an  ion  pair  Eai  = 
0.96  eV  (Eigen  and  De  Maeyer  1956, 1957, 1958; 
Eigen  et  al.  1964),  for  a  long  time  accepted  as  the 
true  value,  should  be  treated  as  the  lowest  possible 
limit  of  this  quantity.  As  the  current  upper  limit  we 
have  Eai = 2E<^ = 1.4  eV  from  the  work  of  Petrenko 
et  al.  (1983).  The  above  discussion  implies  that  ion 
mobility  almost  does  not  depend  on  temperature, 
that  is,  the  ion  motion  is  not  a  thermoactivated 
process.  That  suggestion  was  proven  by  experi¬ 
ment  (see,  e.g.,  Kunst  and  Warman  1983). 

The  determination  of  diarge  carrier  mobilities 
in  ice  is  an  extremely  difficult  task,  since,  because 
of  small  values  of  p;  and  the  "jump"  mechanism  of 
motion,  the  regular  technique  for  measuring  the 
Hall  effect  is  not  applicable  in  this  case.  (The  Hall 
effect  is  discussed  in  the  future  report  on  the 


Surface  of  Ice.)  Hopes  of  using  "saturation  cur¬ 
rents"  for  determining  pi  are  not  justified  because 
of  reasrais  pointed  out  by  Maidique  et  aL  (1971)  in 
dieir  critique  of  Eigen's  work,  and  also  because 
when  high  voltage  is  applied  to  a  system  with  ice  / 
metal  interfaces,  various  physical  phenomena  are 
observed  that  make  the  interpretation  of  results 
more  complicated  (see  Part  II  of  Petrenko  1993a). 

Nonetheless,  several  new  techniques  for  deter¬ 
mining  charge  carrier  mobilities  in  ice  have  been 
developed  and  employed  during  the  last  25  years. 
These  are:  the  technique  that  measures  stationary 
and  transient  infection  ounents  (Eckener  et  al. 
1973,  Petrenko  et  al.  1983,  Petrenko  and  Ryzhkin 
1984b);  the  ice  field  effect  transistor  (Petrenko  and 
Maeno  1987);  and  the  technique  that  measures 
low-frequency  dispersion  (Petrenko  and  Ryzhkin 
1984a,  Zaretskii  et  al.  1988,  Zaretskii  1991).  The 
idea  of  determining  p;  and  nt  from  the  parameters 
of  low-frequency  dispersions  was  presented  by 
Petrenko  (1993a).  Ofoer  techniques  mentioned 
above  are  described  and  discussed  in  Part  n  of  the 
same  review.  The  results  of  determining  charge 
carrier  mobilities  in  ice  are  stunmarized  in  Table  5. 
At  first  glance  only  the  data  on  L-defect  mobilities 
seem  to  be  satisfactory.  At  T= -lO'C,  p4  =  (2  -5)  x 
1(H  m2  V-i  s“i.  This  mobility  drops  exponentially 
with  decreasing  temperature,  the  activation  en¬ 
ergy  being  0.19  -  0.235  eV  (see  Table  6).  There  is 
still  no  reliable  information  about  E)-defect  mobil- 
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Tabk  5.  Mobilities  of  charge  carriers  in  kc. 


T 

H3O* 

W 

(m^A's) 

Oth 

<mWs) 

D 

(h^/Vs) 

L 

(mWs) 

Method  used 

Reference 

-13  .  ..  -36 

(1.1  ±0.1)x  10-7 

— 

— 

— 

Analysis  of 
conductivity  decay 

Kunst  and  Waniian  (1983) 

-10 

2.5  X  10^ 

— 

— 

asandEn 

measureinents 

Wdtz  and  Cole  (1969) 

0...  -20 

5  X  10-7 

— 

— 

— 

Satuiatioa  cunents 

Maidique  et  al.  (1970) 

-123 

— 

54x10-7 

— 

— 

Dofring  with  KOH 

Zarclskii  et  al.  (1988) 

-5...  -40 

10-*-  10-7 

— 

— 

_ 

Proton  injection 

Petrenko  et  al.  (1983) 

0 

2.7  X  10-« 

— 

— 

5xl(H 

Doping  with  HF 

Camplin  et  al.  (1978) 

-10 

— 

— 

— 

l.8xl(H 

Doping  with  HF 

Jaccard  (1959) 

— 

— 

a2x  10-10 

— 

Doping  with  KOH 

Howe  and  Whitworth 
(1989) 

-10 

(2.4  *  1.6)  X  10-7 

— 

— 

— 

Satuiatioo  cunents 

Bullemer  et  al.  (1969) 

-10 

3.5  X  10-* 

— 

— 

— 

HF  doping 

C^aitgilin  and  Glen  (1973) 

-33.1 

9.2  X  10-« 

2.7  X  10-* 

— 

— 

Ice  field  effect  transistor 

Petrenko  and  Maeno 
(1987) 

-145...  -178 

(9±l)xl0-« 

— 

— 

— 

Transient  injection  of 
protons 

Eckener  et  al.  (1973) 

-10 

“  10-* 

— 

— 

l.7xl(H 

Low-fiequency  dispersions 

2!aietskii  et  al.  (1987a) 

-33 

(6.7  ±  0.8)  X  10-7 

— 

— 

__ 

Low-fietpiency  dispersions 

Zaietskii  (1991) 

-20 

6xl0-« 

— 

— 

— 

Proton  injection 

Petrenko  and  Ryzhkin 
(1984c) 

-10 

— 

2.4  X  10-4 

— 

— 

Recombination  injection 

Mtenko  and  Chesnakov 
(1990c) 

-10 

(9±3)xlO-» 

7x10^ 

— 

— 

Analysis  of  data 

Chen  et  al.  (1974) 

ity.  There  are  only  indirect  evidences  that  (13  <  ^4. 
For  excunple,  Devlin  (1992),  citing  spectroscopic 
data,  indicated  that,  in  the  region  at  150  K,  pa  is 
orders  of  magnitude  less  than  P4. 

Data  on  ionic  mobilities  are  distinguished  at 
first  glance  by  their  large  scatter  (from  5  x  10-^ 
down  to  2.5  x  10"*  for  H3O+  and  from  4  x  lO"^ 
down  to  2  X 10-10  for  OH").  However,  if  we  discard 
the  large  values  obtained  by  the  now-discredited 
method  of  saturation  currents  and  the  values  ob¬ 
tained  from  measurements  of  low-frequency  dis¬ 
persions  (which  still  have  low  accuracy),  it  ap¬ 
pears  that  the  positive  ion  mobility  lies  within  the 
limits  of  2.7  x  10"*  and  1.1  x  10"^  m^/V  s — ^notbad 
at  all!  The  authors  of  this  review  prefer  the  data 
obtained  in  transient  injection  experiments 
(Eckener  et  al.  1973),  in  which  a  minimvim  of 
assumptions  was  used.  The  same  H3O+  ion  mobil¬ 
ity  was  foimd  by  Petrenko  and  Maeno  (1987)  using 
field  effect  transistors.  So  it  is  likely  that  pi  =  10*7 
m2  V-i  s"i.  It  seems  well  established  that  both 
H3O  and  OH“  ions  move  in  ice  with  zero  activa¬ 
tion  energy. 

Concluding  this  section  about  point  defects' 
mobilities,  we  can  now  explain  the  at  first  slightly 


surprising  fact  that  the  self-diffusion  coefliaents 
of  oxygen  and  hydrogen  in  ice  are  equal  (Table  2). 
This  is  imexpected  because  the  oxygen  can  diffuse 
only  as  a  part  of  a  water  molecule,  whereas  fite 
hydrogens  can  also  move  by  ion  and  Bjemim 
defect  migration.  Let  us  compare  the  absolute 
values  of  diffusion  coefficients  of  hydrogen  in 
these  two  processes.  If  the  hydrogens  diffuse  in  ice 
by  protonic  defect  transfer,  then  the  diffiision  co¬ 
efficient  Dh  is 

p„-niDl-fn;D; 

"HzO 

where  ni;z  and  Di^  denote  the  concentrations  and 
diffusion  coefficients  of  those  protonic  defects  that 
limit  the  transfer  of  protons  over  macroscopic 
distances.  In  pure  ice  these  defects  are  ions.  There¬ 
fore,  we  can  rewrite  eq  57  in  the  form 

p  _  ^  gs  *  ,5gx 

”  «1^H20  *  "HzO 

where  is  the  static  electric  conductivity  of  ice. 
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For  Dh  to  exceed  the  self-diffusion  coefficient  D*, 
the  following  inequality  must  hold 

S  Dgl— -3°  =  3.9  X  10-*  (n  m)-'  .  (59) 

kgT 

In  eq  59  we  have  used  D*  =  2  x  10-^5  s-t  (T  = 
-10°C)  from  Table  2.  Since  in  pure  ice  at  -10®C 
S  1(H  (n  m)-t  (Petrenko  1993a),  the  diffusion  rate 
of  hydrogen  in  ice  is  faster  when  it  diffuses  as  a 
part  of  a  water  molecule  in  the  form  of  an  intersti¬ 
tial.  That  explains  the  equality  of  the  diffusion 
coefficients  of  hydrogen  and  oxygen.  However,  in 
ice  heavily  doped  with  HF  or  KOH,  would  be 
sufficiently  large  for  hydrogen  transfer  by  pro¬ 
tonic  defects  to  be  faster  fiian  that  by  interstitial 
diffusion  of  H2O  molecules.  There  are  no  experi¬ 
mental  data  on  self-diffusion  coefficients  in  doped 
ice. 

Activation  volume  of  protonic  defects 

The  formation  of  any  protonic  point  defect  re¬ 
sults  in  deformation  of  the  surrounding  lattice.  In 


the  case  of  Bjerrum  defects,  ttie  defective  hydro¬ 
gen  bcmds  are  weaker  than  the  normal  ones,  and 
therefore  they  will  be  longer.  This  is  obvious  in 
particular  for  a  D-defect,  where  ttiere  exists  an 
extremely  strong  repulsicm  betweot  two  protons 
on  one  hydrogen  bond.  Accordingly,  one  can  as¬ 
sume  that  Bjerrum  defects  "expand"  ice,  increas¬ 
ing  its  volume.  If  ice  is  under  pressure  P  then  the 
formation  energy  of  file  defect  Eai  is  increased 
by  the  amount  of  work  against  this  pressure 

£a/fW  =  £a/0J  +  Py„  (60) 

where  Yai  is  the  so-called  activation  volume.  For 
the  defects  fiiat  increase  volume,  Yai  is  positive  and 
their  formation  energy  is  raised  proportionally  to 
the  applied  pressure.  This  results  in  an  exponen¬ 
tially  decreasing  concentration  of  such  defects  with 
increasing  pressure  (Table  7). 

It  is  a  much  more  difficult  problem  to  predict 
even  the  sign  of  Yai  for  ions  in  ice.  It  was  believed 
initially  that,  owing  to  Coloumb  attraction  of  the 
dipole  moments  of  neighboring  water  molecules. 


Table  6.  Data  on  high-frequency  (a^)  and  static  (Os)  conductivities  of  ice. 


T  A*  Of  Esm  Ea 


Type  Of  ice 

{'Q 

(flm)-! 

(Qm)-> 

(eV) 

(eV) 

4  zone  refinings 

-10 

4xl(H 

054 

single  crystals 

— 

2.5xl0-» 

(o...-2(yx:) 

— 

0-0.06 

single  crystals 

-4 ...  -So 

6.4X10-W 

(-10X) 

0.53  ±0.01 

0.70  ±  0.07 

single  crystals 

-10 

— 

10-7 

0.37 

— 

single  crystals 

— 

— 

— 

0.624 

— 

sin^e  crystals 

-10 

1.6x10-5 

1.1  X  10-7 

0.575 

0.61 

pure  aiul  HF-doped 

single  crystals 

0 

— 

2.5  xKH 

— 

— 

single  crystals 

-10 

_ 

(1.1  ±  05)  X 10^ 

0.76  ±0.06 

_ 

single  crystals 

— 

— 

2.3  X  10-* 

— 

0.49 

single  crystak 

-10 

1.8  X 10-5 

— 

0.61 

— 

single  crystals 

-10 

— 

— 

— 

— 

HCl  doping 

S-39 

— 

— 

0.492  ±  0.025 

0 

1(H  mol/L 

HCl  doping 

-39 

— 

— 

0.259  ±0.003 

0.267  ±0.009 

KOH  doping 

0.1  mol/L 

0...  -150 

2xl0-5-10» 

up  to  10-> 

0 

(T<;-70°C) 

0.85 

{T<-70°C) 

(in  solution) 

KOH  doping 

-23...  -203 

-2x10-* 

10-* 

0 

0.4 

0.032  mol/L 

(-23’>C) 

(-23°C) 

(T>-100^) 

(T>-100^) 

(in  solution) 
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ions  will  contract  the  lattice,  i.e.,  Yai  <  0.  Experi¬ 
mental  data  shown  in  Table  7  refuted  this  opinion, 
and  the  activation  volumes  of  both  riaO  and  OH~ 
iotts  turned  out  to  be  positive.  We  can  try  to 
explain  this  by  the  fact  that  dipole  moments  of  all 
neighboring  molecules  in  the  ice  structure  cannot 
be  oriented  towards  the  ion.  This  would  result  in 
a  violation  of  the  Bemal-Fbwler  rules.  Also,  if  a 
dipole  moment  is  directed  away  from  the  ion,  it 
will  be  repelled  rather  than  attracted,  increasing 
the  voliune  around  the  def^.  Thus,  such  simple 
considerations  are  not  sufficient  to  estimate  Yai  of 
ions,  and  ab  initio  calculations  taking  into  accotmt 
the  permitted  configurations  of  water  molecule 
orientations  in  ice  are  required. 

For  a  defect  to  jump  from  one  molecule  to  the 
neighboring  one  (or  from  one  bond  to  another 
one),  an  additional  change  in  voliime  Ymi  may  be 
required.  This  is  rather  obvious,  since  the  energy 
of  the  transitional  state  of  a  defect  is  of  course  not 
lower  than  the  energy  of  the  equilibrium  state. 
Therefore,  when  coirsidering  a  defect  transfer  im- 


derpressure,  we  have  to  add  an  additirmal  term  to 
the  activation  energy  of  their  motion  Emi 

^Bu(P)  —  Pmi(0)  +  Tfc(u  P .  (61) 

In  other  words,  the  amplitude  of  the  potential 
relief  shown  in  Figure  1  varies  with  pressure. 
During  ion  transfer,  whidr,  according  to  experi¬ 
mental  data,  occurs  through  a  non-activatioiral 
hop  of  protcms  along  a  bond,  the  crucial  factor  may 
be  an  exponential  increase  in  the  probability  of 
such  a  transition  attributable  to  a  decrease  in  the 
length  of  the  hydrogen  bond  owing  to  the  applied 
pressure.  This  would  appear  as  a  negative  activa¬ 
tion  volume  Ymi- 

In  ffie  course  of  studies  on  the  dependence  of 
partial  conductivities  of  defects  Oi  on  pressure,  it 
was  found  that 

oi=Nm«i‘>‘ 

^ai(Q)  ^mi(Q)]  ^Yai  Yiiii]\  /^2) 


Table6(cont'd). 

Activation 


(eV) 

EoB 

(eV) 

Activation 
energy  of 
ion  mobility 
(eV) 

energy  of 

D-  and  Irdefeci 
mobility 
(eV) 

td 

<s) 

E, 

(eV) 

Reference 

1.08 

_ 

_ 

WOiz  and  Cole  (1969) 

l.l 

— 

— 

— 

5.18  X  10-S 

0.598 

Maidiqije  et  al.  (1970) 

(-10"C) 

1.4 

— 

— 

— 

— 

— 

Pebenko  et  al.  (1983) 

__ 

_ 

_ 

BuUemer  and  RkM  (1966) 

— 

— 

— 

— 

5x10-5 

— 

Camplin  et  al.  (1978) 

1.22 

0.68 

0 

0.235 

2x  10-’ 

0.575 

Jaccard  (1959) 

(0°Q 

— 

— 

— 

— 

1.8x10-5 

0.59 

von  Hippeletal.  (1971) 

Maidique  et  al.  (1971) 

— 

— 

— 

— 

— 

— 

Bullemer  et  aL  (1%9) 

0.98 

— 

— 

— 

Maeno  (1973) 

_ 

— 

— 

— 

4.5  X  10-5 

0.64 

Camplin  and  Glen  (1973) 

_ 

— 

— 

— 

10-4 

0.57 

Zaietskii  et  al.  (1989) 

_ 

— 

0 

— 

— 

0.567  ±0.00 

Takei  and  Maeno  (1984) 

90.293  ±  0.034 

_ _ 

0.79 

0 

0.19 

— 

— 

Take!  and  Maeno  (1987) 

_ 

_ 

0 

— 

2x10-5 

0 

Zaretskii  et  al.  (1988) 

(7<-70‘’C) 

_ 

Howe  and  Whitworth 

0 

(1989) 
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Table  7.  Activation  volumes  of  protonic  point  defects  in  ice.  The  data  are  given  in  the  units  of  V-.~i 
=  3.3  X 10-23  cm3. 


/ce  type 


(concentrations 
in  mole/L) 

yet 

(H}0*  ion) 

y<a 

(OH-  ion) 

yaJ 

(D-defect) 

y»4 

(Lr  defect) 

yel.2  +  ymU 

YaJ,4  +  ymJ,4 

Reference 

pure  polycrystais 

— 

— 

— 

— 

-{0.55  ±0.15) 

0.145  ±0.015 

Chen  eial.  (1965) 

pure  single  crystals 

-(0.18  ±0.01) 

0.13  ±0.01 

Taubenberger  et  al. 
(1973) 

-  5  X  10-4  HF 

0.091  ±  0.024 

-(0.16  ±  0.024) 

Taubenberger  et  al. 
(1973) 

-3x  10-5HF 

— 

0 

-(0.36  ±0.05) 

Taubenberger  et  al. 
(1973) 

pure  single  crystals 

— 

— 

— 

— 

-{3  ±0.9) 

•  0.15 
(-25"C) 

Hubmann  (1978) 

1(HNH3 

— 

— 

— 

— 

--0.7 

-0.27 

(-TS’O 

Hubmann  (1978) 

3  X  10-«  HF 

— 

— 

— 

--0.55 

-0.15 

(-25°0 

Hidimann  (1978) 

single  crystals 

KM  Ha 

3.3  ±0.12 

— 

1.33  ±0.3 

— 

— 

Evtushenko  and 
Petrenko  (1991) 

single  crystals 
KHHF 

3.9  ±0.36 

— 

— 

1.0  ±0.39 

— 

Evtushenko  and 
Petrenko  (1991) 

single  crystals 

1(H  NH3 

" 

0.82  ±0.39 

0.7  ±0.3 

■ 

■ 

Evtushenko  and 
Petrenko  (1991) 

We  can  determine  only  the  "total"  activation  vol¬ 
ume  [Yai  +  Ymil-  Such  measurements  were  origi¬ 
nally  performed  by  Chen  et  al.  (1965)  on  pure 
polycrystalline  ice.  They  fotmd  that  the  static  con¬ 
ductivity  Os  increases  and  the  high-frequency  con¬ 
ductivity  (T„  decreases  with  increasing  pressure. 
Similar  measurements  on  monocrystalline  ice  were 
performed  later  by  Taubenberger  et  al.  (1973)  and 
Hubmann  (1978).  Taubenberger  et  al.  used  pure 
and  HF  doped  ice  and  Hubmann  used  pure  and 
NH3  and  HF  doped  ice.  The  idea  of  using  the 
doped  ice  was  to  produce  a  known  initial  concen¬ 
tration  of  various  defects  to  separate  out  Yai  md 
Ymi.  However,  the  interpretation  of  the  results 
obtained  from  the  doped  ice  still  remains  compli¬ 
cated  and  ambiguous,  since  pressure  can  affect  the 
dissociation  energies  of  protonic  defects  from  im¬ 
purities.  Table  7  contains  values  of  (Yai + Jm)  found 
in  these  works.  Here,  we  assume  that  and 
are  proportional  to  the  partial  conductivity  of 
the  majority  charge  carrier  and  that  Os  is  propor¬ 
tional  to  the  partial  conductivity  of  the  minority 
charge  carrier,  in  conformity  wi  A  eq  49-51 .  In  the 
analysis  of  data  presented  in  Table  7,  we  have  to 
take  into  accoimt  that  in  the  case  of  doped  ice  Eai  is 
related,  obviously,  to  the  energy  of  separation  of  a 
defect  from  an  impurity. 


Evtushenko  and  Petrenko  (1991)  used  a  quite 
different  approach  to  determine  activation  vol¬ 
umes  of  protonic  defects  in  ice.  Since  ttie  formation 
energy  of  defects  £»  depends  on  the  applied  pres¬ 
sure  P  according  to  eq  57,  in  %e  presence  of  a 
pressure  gradient  there  is  a  force  F,  acting  upon  the 
defects 

Fi  =  -grad(£a, )  =  -Yai  grad(P) .  (63) 

In  their  work,  Evtushenko  and  Petrenko  intro¬ 
duced  a  pressure  gradient  by  bending  thin  ice  sam¬ 
ples  and  then  measured  the  electric  field  sbength 
E  compensating  for  the  effect  of  the  force  Fj 

eiE  =  -yafgradfP).  (64) 

By  varying  temperature  and  doping,  they  deter¬ 
mined  activation  volumes  of  ^  four  t^es  of 
protonic  defects  in  ice  (see  Table  7).  The  advantage 
of  this  method  is  that  it  enables  us  to  separate  out 
£ai  (i.e.,  the  additional  term  to  the  specific  voliune 
of  a  defect)  from  the  sum  £ai + Em/.  The  results  pre¬ 
sented  in  the  table  are  expressed  in  terms  of  the 
volume  occupied  by  a  water  molecule  in  ice:  3.3  x 
10-23  cm3.  According  to  the  results  of  Evtushenko 
and  Petrenko,  all  protonic  defects  (except  HaO'*' 
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i(m)  have  approximately  doubled  volume  com- 
paied  to  die  molecular  volume. 

Spectroscopic  study  of  protonic  dcf^ 
migration — ^D-H  exchange 

An  original  and  informative  optical  method  of 
investigating  protonic  defects  in  ice  was  devel¬ 
oped  by  Devlin  (Ritzhaupt  and  Devlin  1977; 
Ritzhaupt  et  al.  1978,  1980;  Devlin  1990).  The 
method  is  based  on  the  optical  detection  of  hydro¬ 
gen  isotope  exchange  (H  and  D)  between  water 
molecules  in  ice.  Since  such  an  exchange  can  occur 
only  through  protonic  defect  transfer,  the  rate  of 
isotope  exdumge  enables  tis  to  estimate  the  di^- 
sion  of  ions  and  Bjemim  defects  in  ice.  Figure  8 


illustrates  the  essence  of  this  method.  If  we  place  a 
molecule  in  H20ke,  then,  because  of  the  trans¬ 
fer  of  ions  displacing  protcms  and  deuterons  akmg 
hydrogen  bonds,  pairs  of  HOD  molecules  will  be 
formed  on  adjacent  sites.  We  shall  denote  these 
pairs  as  (HOD)2.  Hence,  the  rate  of  the  reaction 

D2O  +  H2O  ->  (H0D)2  (65) 

will  be  determined  by  die  number  of  intermolecu- 
lar  jumps  by  ions,  which  is  proporticmal  to  d\e 
ionic  conductivity 

For  two  deuterons  to  move  a  large  distance 
apart,  it  is  necessary  to  transfer  a  deuteron  from 
one  bond  to  another.  This  can  take  place  (xily 


Exan^e  of  ion 
and  L-  defecf 


I 


H 


Turn  step 

L-defect  migration  which 
separates  the  D  by  another 
oxygen. 


Figure  8.  Point  d^ts  in  an  ice  lattice,  showing  the  proton  hop  step  and  the  L-defect  turn  step 
of  proton  transport  and  isotopic  exchange. 
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through  Bjemun  defect  migration.  Therefore,  the 
rate  of  the  reaction 

(H0D)2  ->  HOD  +  HOD  (66) 

will  be  determined  by  the  number  of  elementary 
jumps  of  Bjemun  defects  per  second,  which  is 
proportional  to  pB  <%• 

The  reactions  (eq  65-66)  can  be  observed  be¬ 
cause  the  infrared  absorption  spectra  of  D2O, 
(H0D)2  and  isolated  HOD  are  different,  as  shown 
in  Figure  9.  Devlin  and  co-authors  prepared  H2O 
ice  samples  with  a  small  addition  of  D2O  by  vapor 
deposition  in  a  vacuum  at  125  K.  At  such  a  low 
temperature,  the  protonic  defect  transfer  is  negli¬ 
gibly  small  and  the  ice  obtained  has  the  cubic 
structure  Ic-  Varying  the  temperature  in  5-K  steps 
and  observing  the  kinetics  of  the  changes  in  the  IR 
spectra,  E)evlin  and  co-authors  (Collier  et  al.  1984) 
found  that  in  their  samples  in  the  temperature 
range  from  130  to  150  K 


Figure  9.  Infrared  spectra  for  isotopomers  ofHzO 
decoupled  in  H2Oiceat90  K:  (a)  D2O;  (b)  (HOD)2; 
(c)  HOD  (after  Devlin  1992). 


Eai*0.41eV 

and  (67) 

EoB»0.52eV 

which,  although  close  to  the  values  determined 
from  electrical  conductivity  (see  Table  6),  are  still 
considerably  lower  than  the  conventional  values 
at  higher  temperatures.  This  discrepancy  may  be 
attributed  to  a  number  of  factors,  including  extrin¬ 
sic  defect  activity,  the  cubic  natiue  of  the  ice  used, 
or  the  predicted  behavior  ignoring  shallow  proton 
trapping,*  or  all  three. 

Probably  the  major  significance  of  the  work  of 
Devlin  and  colleagues  is  that  they  confirmed  at  the 


*The  shallow  proton  traps  were  identified  with  the  nega¬ 
tively  charged  L-defects  (Kunst  and  Warman  1983).  An  asso¬ 
ciation  energy  of  a  proton  on  such  a  trap  was  found  as  «0.43  eV 
by  Woodbridge  and  Devlin  (1988)  and  >0.19  eV  by  Kunst  and 
Warman  (1983). 


Figure  10.  Oxygen  atom  positions  in  ice  Ih. 
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molecular  level  the  two-step  character  of  protcm 
trai\sfer  in  ice,  and  consequently  the  necessity  for 
such  a  transfer  to  be  performed  by  two  groups  of 
very  different  defects:  ions  and  Bjerrum  defects.  In 
other  words,  the  work  of  this  group  provides  very 
strong  support  for  the  existence  of  both  t)rpes  of 
protonic  defects  in  ice.  Readers  interested  in  this 
question  can  find  details  in  the  review  by  Devlin 
(1990).  We  shall  also  return  to  this  topic  in  the 
report  Optical  Properties  of  Ice. 

Atomic  structure  of  protonic  defects 

The  atomic  structure  of  protonic  defects  de¬ 
picted  in  Figure  4b  is  purely  schematic.  The  intro¬ 
duction  of  a  defect  changes  considerably  the  forces 
acting  between  adjacent  atoms,  and  this  leads  to 
lattice  "relaxation"  aroimd  the  defect.  One  piece  of 
evidence  for  such  relaxation  is  the  increased  spe¬ 
cific  volume  of  defects  as  compared  to  the  molecu¬ 
lar  volume  (see  the  Activation  Volume  of  Protonic 
D^ects  section).  The  necessity  for  atomic  rear¬ 
rangement  of  the  ice  lattice  around  protonic  de¬ 
fects  can  be  demonstrated  by  estimation  of  the 
electrostatic  energy  of  interaction  between  two 
protons  on  the  same  hydrogen  bond  in  a  D-defect. 
If  we  do  not  take  relaxation  and  lattice  rearrange¬ 
ment  into  account,  then  the  two  protons  will  be  at 
a  distance  r  of  about  0.7  A  from  each  other.  Since  at 
such  smaU  distances  we  have  to  use  c  =  1,  the 
electrostatic  energy  of  such  a  D-defect  would  be 

AE«-e^=19eV.  (68) 

\mreQ 

This  is  extremely  large  compared  to  the  observed 
formation  energy  of  0.68  eV  of  a  peiir  of  L-  and  D- 
defects. 

Currently,  there  are  two  models  in  which  the 
energy  of  protonic  defects  is  lowered  by  rew- 
rangement  of  atoms  around  a  defect.  In  the  first 
model,  after  a  protoiuc  defect  is  formed  in  the 
simplified  fashion  shown  in  Figure  4b,  atorrrs  and 
molecules  are  2dlowed  to  shift  sbghtly  from  their 
initial  positions.  The  new  equilibrium  stable  con¬ 
figuration  is  the  one  with  the  lowest  energy.  This 
approach  to  the  calculation  of  energy  and  struc¬ 
ture  of  Bjerrum  defects  in  ice  was  used  by  Plummer 
(1987, 1W2).  She  used  ab  initio  computer  simula¬ 
tions,  requiring  very  large  amounts  of  computer 
time,  and  because  of  this  only  clusters  comprising 
eight  water  molecules  with  a  Bjerrum  defect  in  the 
middle  could  be  simulated.  Two  ways  of  choosing 
such  clusters  with  "eclipsed"  and  "staggered" 
structure  are  shown  in  Figure  10.  The  new  inter¬ 
atomic  distances  and  energies  AE  for  the  formation 


Table  8.  Calculated  parameters  of  Bjerrum 
defects  (after  Plummer  1992). 

Parameter  Eclipsed  Stagger^ 


D-defect 


0-0  distance  (A) 

3.30 

2.72 

H-H  distance  (A) 

1.78 

12i 

Creation  energy  (eV) 

0.417 

0.278 

L-d^ect 

0-0  distance  (A) 

3.18 

3.16 

Creation  energy  (eV) 

0.317 

0.343 

of  L-  and  D-defects  in  such  clusters  are  summa¬ 
rized  in  Table  8. 

Since  in  a  real  ice  lattice  Bjerrum  defects  can  be 
located  at  hydrogen  bonds  Irath  parallel  (eclipsed 
structure)  and  inclined  (staggered  struchire)  to  the 
hexagonal  c-axis,  we  should  take  some  average 
between  the  eclipsed  and  staggered  values  for 
comparison  with  experimental  values,  taking  into 
account  that  there  are  three  times  more  bonds 
inclined  to  the  c-axis  than  are  parallel  to  it.  The 
results  of  this  averaging  are  given  in  Table  9,  along 
with  values  for  the  ideal  ice  structure  Ih,  but  we 
must  remember  that  the  defects  have  been  mod¬ 
eled  not  in  bulk  ice  but  in  a  cluster  of  only  eight 
molecules. 

The  first  thing  that  attracts  our  attention  is  fiiat 
the  lengfii  of  the  bond  carrying  the  defect  has 
incrcas^.  This  increase  is  larger  for  L-defects, 
which  is  in  qualitative  agreement  with  the  larger 
activation  volume  for  L-defects  foimd  experimen¬ 
tally  (Table  7).  Plummer  discovered  a  large  in¬ 
crease  in  the  H-H  distance  in  a  D-defect.  One  of  the 
protons  actually  moves  off  the  hydrogen  bond, 
taking  an  interstitial  position,  and  the  remaining 
hydrogen  bond  becomes  bent  in  the  opposite  di¬ 
rection  by  approximately  20®. 

It  is  quite  remarkable  that,  in  spite  of  the  small 
size  of  simtdated  clusters,  the  energy  of  pairs  of  L- 


Table  9.  Average  parameters  of  Bjerrum  defects 
in  ice. 


Parameter 

Defect 

Perfect 

structure 

D 

L 

0-0  distance  (A) 

2.865 

3.17 

2.76 

H-H  distance  (A) 

2.12 

- 

0.76 

Creation  energy  (eV) 

0.312 

0.336 

- 
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a.  Ionic  de^t  in  one^imensional  chain  (^hfdrogen  bonds. 


Figure  11.  Soliton  model  of  OUT  ion. 


and  EMefects  according  to  Plummer  is  about  0.65 
eV,  which  is  very  close  to  the  value  0.68  eV  found 
from  experiments. 

We  can  make  one  more  interesting  conclusion 
from  Pluixuner's  calculations  by  observing  that 
energies  of  defects  in  "eclipsed"  and  "staggered" 
configurations  are  different.  Because  in  diffusion 
and  electrical  conductivity  processes,  Bjerrum 
defects  must  flip  from  one  configuration  to  an¬ 


other,  the  difference  in  corresponding  energies 
may  contribute  to  the  observed  activation  energy 
of  defect  migration.  For  L-defects  this  difierence  is 
ordy0.026eV  (i.e.,  on  the  order  of  IbaT),  whereas  for 
I>defects  it  is  equal  to  0.14  eV.  Is  this  not  the  reason 
for  lower  mobility  of  l>defects  compared  to  L- 
defects?  Calculations  of  the  type  poformed  by 
Pltmuner  are  quite  complicated  and  have  not  yet 
been  made  for  ions. 


■WPo 


a.  Defect  in  one-dimensional  chain. 


Figure  12.  Soliton  model  of  Bjerrum  defect. 
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The  second  type  of  model  of  dte  atomic  struc¬ 
ture  of  protonic  defects  is  the  so-called  soliton 
model,  in  which  the  def^  is  spatially  delocalized 
over  several  lattice  sites.  Both  ekctrostatic  and 
elastic  energies  of  the  defect  are  lowered  by  this. 
Origiiudly  suggested  for  the  description  of  ion 
transfer  along  hydrogen  bond  chains  by 
Antcmchenko  et  aL  (1983)  and  for  E^rrum  defects 
by  Sergienko  (1986),  this  model  has  been  furdier 
developed  in  many  papers;  it  has  numerous  advo¬ 
cates  and  no  less  numerous  opponents.  A  detailed 
theoretical  description  of  this  model  can  be  found 
in  a  recent  book  by  Davydov  (1991).  A  simplified 
diagram  of  an  OH~  "soliton  ion"  in  a  linear  chain 
of  hydrogen  bonds  is  shown  in  Figure  11.  Notice 
that  all  protons  to  the  left  of  the  center  are  dis¬ 
placed  to  the  left  end  of  their  bonds  and  all  protons 
on  the  right  are  displaced  to  the  right.  During  the 
migration  of  such  an  ion  to  the  right,  the  protons 
on  the  Ixxids  are  transferred  to  dte  left  For  a  soliton 
to  arise  and  remain  stable,  at  least  two  conditioits 
must  be  satisfied:  1)  the  protons  must  behave 
cooperatively,  i.e.,  they  must  strongly  interact  with 
each  other,  and  2)  a  stroirg  nonlinear  interaction  of 
the  proton  and  oxygen  sub-lattices  is  required. 

The  soliton  simulating  an  L-defect  (the  migra¬ 
tion  of  the  defect  from  left  to  right  transfers  the 
protons  from  right-side  positions  on  the  bonds  to 
left-side  positions)  is  shown  in  Figure  12.  Using 
the  same  schematic,  the  reader  can  easily  construct 
an  HaO*’  ion  and  a  D-defect.  Since  the  transfer  of 
soliton  defects  results  in  exactly  the  same  changes 
in  the  positions  of  protcms  on  the  bonds  as  in  fire 
case  of  transfer  of  their  "point"  analogs,  the  formal 
description  of  their  effe^  in  the  ice,  effective 
charges,  contributions  to  D ,  etc.,  will  be  exactly  the 
same. 

In  spite  of  a  certain  elegance  of  the  soliton 
model,  serious  arguments  have  appeared  against 
its  applicability  to  ice.  The  first  eirgument  pre¬ 
sent^  by  Na^  (1992)  is  that  solitons,  requiring 
for  their  motion  a  coordinated  displacement  of 
protorrs  2uranged  into  some  ordered  positions, 
will  just  get  "stuck"  in  a  disordered  ice  protonic 
structure.  A  second  objecticm  is  that  the  soliton  is 
defined  and  has  to  move  along  a  particular  track, 
whereas  the  bonds  in  ice  form  a  branched  net¬ 
work.  The  relaxed  state  of  a  static  defect  may  in 
reality  be  somewhat  delocalized,  but  the  theory  of 
solitons  is  essentially  dyruunic. 

A  third  objection  to  the  soliton  model  in  or¬ 
dered  water  molecule  chains  arises  from  precise 
ab  initio  calculations  performed  by  Godzik  (1990) 
and  Scheiner  (1991).  These  calculations  show  that 


in  ordered  chains  of  H2O  mc^ecuks  the  coopera¬ 
tive  motitm  of  several  protons  involves  a  sigrtifi- 
cantly  higher  energy  barrier  fium  firat  for  succes¬ 
sive  jumps  of  one  protcm  after  another. 

We  cannot  leave  the  structure  of  protortic  de¬ 
fects  without  mentioning  the  possibility  that  these 
defects  could  be  coupled  wiA  vacancies  (Kroger 
1974).  For  example  a  I>-defect  at  a  vacarwy  is 
merely  a  vacant  site  with  three  rafirer  than  the 
norm^  two  hydrogen  Ixnids  pointing  into  it,  and 
its  D  character  will  move  with  it  as  it  moves 
through  the  lattice  (see  Fig.  17).  As  vacancies  may 
be  highly  mobile,  we  must  allow  for  the  possibility 
that  the  majority  of  Bjerrum  defects  could  exist 
and  move  in  this  way.  There  is,  however,  no  ex¬ 
perimental  evidence  on  this  issue. 


IMPURITIES 

Althou^  liquid  water  is  a  very  good  solvent  for 
many  substances,  its  solid  phase — ^ice — dissolves 
almost  nothing.  Most  of  fiw  impurities  dissolved 
in  water  either  get  rejected  as  fiw  water  freezes  or 
precipitate  as  a  second  phase  within  the  ice.  There 
are  some  exceptiorrs,  such  as  a  few  adds  (HF,  HCl), 
ammonia  (NH3),  some  alkalis  (KOH,  NaOH)  and 
their  derivatives  like  NH4F.  These  substances  can 
be  incorporated  properly  into  the  ice  structure 
where  firey  drastically  change  fiie  protonic  carrier 
concentrations  in  fire  ice.  To  understand  how  these 
impurities  afreet  the  electrical  properties  of  ice,  let 
us  emrsider  the  probable  patterns  of  their  incorpo¬ 
ration  into  the  ice  lattice  (see  Fig.  13-15).  Note  ^t 
alfirough  there  now  exist  a  number  of  experimen¬ 
tal  techniques  for  investigating  crystal  structures 
on  an  atorrtic  scale  (high  resolution  transmission 
electron  microscopy,  scanning  tunneling  micros¬ 
copy,  atomic  force  microscopy,  etc.),  such  tech¬ 
niques  have  not  yet  been  successfully  applied  to 
ice.  Patterns  such  as  those  represented  in  Figures 
13-15  turn  out  to  be  just  hypofiieses  based  upon 
ionic  radii,  energies  and  lengfiis  of  atomic  bonds, 
etc.  The  ionic  radius  of  F~,  for  example,  is  rather 
dose  to  that  of  O^-,  and  it  is  natural  to  assume  that 
fluorine  substitutes  for  oxygen  successfully  in  an 
ice  lattice.  As  can  be  seen  from  Figure  13,  since  the 
HF  molecule  can  provide  only  one  proton  for  the 
four  nearest  hydrogen  bonds,  one  bond  appears  to 
be  lacking  any  protons  at  all;  i.e.,  substitutiem  of 
the  water  mole^e  by  the  HF  molecule  results  in 
the  appearance  of  an  L-defed  associated  with  the 
impimty  molecule.  The  binding  energy  of  such  an 
L-defed  with  the  HF  molecule  is  not  large  (^  10-^ 
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eV  (Jaccaid  1959,  Camplin  et  al.  1978]).  Since 
number  of  possible  arrangements  of  the  L>defect 
in  an  ice  lattice  is  enormous,  practically  all  L- 
defects  leave  their  pamtt  molecules. 

Moreover,  it  turns  out  that  the  single  proton  in 
the  HF  molecule  is  not  bound  firmly  and  can  jump 
along  a  hydrogen  b<md  to  anofirer  water  molecule, 
turning  drat  water  molecule  into  an  H3O  ion.  It  is 
important  to  emphasize  the  dwmistry  of  HF.  To  a 
physicist  the  rratural  thing  would  be  for  an  extra 


proton  to  attach  to  HF,  forming  an  OH~  icm,  but 
chemically  HF  is  a  weak  acid  in  water  and  prefers 
to  lose  a  proton. 

At  the  same  time  as  HF  doping  of  ice  increases 
the  L-defect  and  HaO  icm  concentratioirs,  it  will 
also  lower  the  concortrations  of  D-defects  and 
OH-  ions.  This  happetts  because  the  equilibrium 
ccmcentraticHis  of  D^efects  and  OH-  ions  in  pure 
ice  are  determiired  by  a  balance  betweor  fire  rate  of 
fireir  thermoactivat^  production  and  the  rate  of 


b.  Same  after  HsO*  ion  and  L-defecl  were  released  by  means  of  sequential  hops  of  protons 
(along  the  hydrogen  bonds  to  release  H3O*  ion  and  between  the  bonds  to  release  L-defect). 

Figure  13.  Hypothesized  pattern  of  incorporation  ofHF  impurity. 


recombination:  D-  with  L-defects  and  H3O  with 
OH~  ions.  HF  doping  increases  L-defect  and  HaO^ 
ion  concentrations,  so  that  the  rate  of  recombina¬ 
tion  increases,  which  results  in  decreasing  I>de- 
fect  and  OH-  ion  concentraticms.  This  can  be  stated 
in  terms  of  the  law  of  mass  action,  whidi  is  that  the 
products  of  concentrations  ni  x  n2  and  na  x  114  must 
be  constant  in  thermal  equilibrium,  irrespective  of 
doping 

«l  X  «2  =  I  («H2o)^xp  ,  nij2  *  nH20 


a.  NH3  molecule  incorporated  into  ice  structure. 


b.  Same  after  release  ofOH~  ions  and  D-defect  from  the  NH3  molecule. 


Figure  14.  Hypothesized  pattern  of  incorporation  ofNH3  impurity. 


M3  X  «4 . 


(«H2ofexpj-i 


«3,4  «  «H20  ■ 


Another  molecule  that  can  successfully  be  in¬ 
corporated  into  the  ice  lattice,  substituting  for 
water  molecules,  is  ammonia,  NH3  (Hg.  14).  Since 
NH3  donates  duee  protons  (instead  of  two  as  in 
the  case  of  H2O)  for  four  hydrogen  bonds,  we 
observe  two  protons  at  (me  Ixmd  at  the  same  time, 
i.e.,  a  D-defect.  As  l.  consequence  of  the 
diermoactivation  process,  dtis  defect  can  be  re¬ 
leased  and  take  p^  in  the  conductivity  piocess. 
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Besides,  the  NH3  molecule  is  a  base  and  can  accept 
a  fourth  proton,  forming  an  OH“  ion  at  a  neighbor¬ 
ing  water  molecule.  Activation  energies  for  the 
release  of  D-defects  and  OH~  ions  from  the  ammo¬ 
nia  molecule  in  ice  are  comparatively  large,  so 
that,  although  ammonia  doping  of  ice  increases  D- 
defect  and  OH"  ion  concentrations  (decreasing  L- 
defect  and  HaO^  ion  concentrations  in  accordance 
with  eq 69-70),  its  influence  in  terms  of  conductiv¬ 


ity  is  not  as  efficiort  as  that  of  HF.  It  is  also 
important  that  NH3  doping  increases  the  concoi- 
tration  of  less  mobile  carriers  (D-defects  and  OH~ 
ions);  see  Table  6. 

Lastly,  the  third,  very  important  group  of  im¬ 
purities  that  consider2ibly  affect  the  conductivity 
of  ice  are  the  alkali  hydroxides,  among  which  the 
most  sigiuficantisKOH.  The  mostplausiblemecha- 
rusm  of  incorporation  of  KOH  into  the  ice  struc- 


a.  How  KOH  dissolved  in  ice  is  thought  to  introduce  defects.  The  K*  ion  is  interstitial,  the 
OH~  ionsubstitutes for  awatermolecukleamngonebo^  without  any  proton— an  L-d^ect. 


b.  Same  after  L-defect  and  OH~  ion  have  moved  axoay  from  the  ff*"  ion. 
Figure  15.  Hypothesized  pattern  of  incorporation  of  KOH  impurity. 


Figure  16.  Absorption  spectra  cf  solvated  electrons  in 
water  and  ice  (after  Nilsson  et  al.  1972), 


ture  is  shown  in  Figure  15.  The  potassium  atcm  is 
accommodated  in  interstices,  while  die  hydroxyl 
group  OH  is  built  into  the  lattice.  As  can  be  seen 
from  Figure  15,  this  results  in  the  appearance  of  L- 
defects  and  OH-  ions  linked  to  KOH  impurity. 
Thermoactivation  processes  allow  them  to  be  re¬ 
leased  into  the  ice  bulk,  increasing  L-defect  and' 
OH-  ion  concentrations  and  suppressing  D-de- 
fects  and  H3O+  ions.  Because  of  the  extremely 
small  binding  energies  of  L-defects  and  OH-  ions 
with  ions  in  ice,  these  impiuities  become  effi¬ 
cient  donors  of  L-defects  and  OH-  i<ms  even  at 
very  low  temperatures  (2^retskii  et  al.  1988). 

The  subject  of  the  concentrations  of  charge  car¬ 
riers  in  doped  ice  has  been  treated  many  times  in 
the  scientific  literature  (see  Jaccard  1959,  Krdger 
1974,  Camplin  et  aL  1978),  and  is  described  more 
fully  in  the  report  on  the  Electrical  Properties  of  Ice 
(Petrenko  1993a).  For  the  fullest  detail,  based  on 
the  law  of  mass  action,  the  reader  is  referred  to 
Kroger  (1974). 

At  present  there  are  himdreds  of  papers  avail¬ 
able  in  the  literature  in  which  the  electrical  proper¬ 
ties  of  doped  ice  were  studied.  The  results  of  these 
studies  were  discussed  in  detail  in  the  review  by 
Petrenko  (1993a)  and  are  presented  in  part  in 
Tables  4-6.  Let  us  refer  to  die  following  works  as 


examples:  Takei  and  Maeno  (1984, 1987)  cm  HCl 
doped  ice;  Gross  (1975)  and  Gross  et  aL  (1978)  on 
HCl,  HF,  NaCl,  KF,  NH4F,  ISfi^OH,  NH4CI, 
NH4HCC>3  d<^)ed  ice;  KOH  doped  ice  was  stuclied 
by  Zaretskii  et  aL  (1988)  and  by  Howe  and 
Whitworth  (1989). 

On  the  whole,  ffie  effect  of  impurities  dissolved 
in  ice  on  its  electrical  properties  is  rather  well 
explained  on  the  basis  of  the  currently  accepted 
model  of  electrical  conductivity  of  ice,  provided 
that  the  concmtrations  of  all  four  types  of  charge 
carriers  are  calculated  following  the  methcxl  sug¬ 
gested  by  Kroger  (1974). 

The  impurity  diffusion  coefficients  in  ice  were 
measured  in  a  number  of  works.  Haltenorth  and 
Klinger  (1%9, 1977)  investigated  diffusion  of  HF 
in  mcmocrystalline  ice  parallel  and  perpendicular 
to  the  c-axis.  They  found  ffiat  at  T  =  -10°C  foe 
diffiision  coefficient  of  HF  in  ice  is  1.08  x  10*it  ii|2 
sri  and  obeys  the  Arrheiuus  law  with  an  activation 
energy  of  0.2  eV.  The  value  of  the  diffiision  coeffi¬ 
cient  along  the  c-axis  was  20%  lower  than  the 
corresponding  value  normal  to  it.  In  polycrystal¬ 
line  ice  foe  diffiisicm  coefficient  was  25%  higher 
ffian  in  mcmocrystals  owing  to  grain  boundary 
diffiision.  It  should  be  noted  that  file  diffusion 
coefficient  for  HF  is  much  larger  than  the  self- 
diffiision  coefficients  of  H2O  given  in  Table  2, 
which  suggests  that  HF  molecules  move  through 
interstitial  positions  more  easily  than  those  of 
H2O. 

Diffusion  of  the  neutral  gases  He  and  Ne  was 
studied  by  Haas  et  al.  (1973).  Diffiision  coefficients 
of  these  gases  in  ice  appeared  to  be  extremely  high, 
indicating  a  "transparency"  of  ice  for  these  gases. 
It  is  interesting  that  the  diffusion  coefficients  along 
the  c-axis  (D|  ic)  were  six  times  (He)  and  eight 
times  (Ne)  hi^ier  fiian  in  the  perpendicular  direc- 
ticm.  In  file  temperature  range  173-253  K  for  He 
and  213-263  K  for  Ne,  Haas  and  co-aufiiors  found 
that  the  foUowing  relationships  hold 

D|  ic(He)  =  3.4  x  lOr^  exp  f-  012  gV\  ^  g-i  (yj) 

\  ksT  I 

D^He)  =  1.1  X  10-7  exp  (-  irf  s-i  (72) 

\  kbT  / 

D|  ic(Ne)  =  1.2  X  l(r6  exp  (-  Q.25_eV  j  g-i  (73) 

\  kbT  I 

ELECTRONIC  DEFECTS 

In  this  section  we  will  consider  point  defects 
that  arise  from  an  excitation  of  the  electronic  sub¬ 
system  of  ice:  hydrated  electrons,  radicals  and 
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Figure  17.  Complex  (mcmcy  +  D-defect)  formation  in  ice  (after  de  Haas  et  al.  1983). 
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color  centers.  Since  ice  has  a  very  wide  forbidden 
band  with  Eg  »  8.5  eV  (Ryzhkin  1992),  the  above 
excitatioirs  usually  arise  when  ice  is  either  illumi¬ 
nated  by  hard  ultraviolet  radiation  or  is  excited  by 
beams  of  charged  particles,  y-rays  and  X-rays. 

Conducting  electrons  in  ice  have  been  studied 
m  many  works  (e.g.,  Verbeme  et  aL  1978;  Warman 
et  al.  1980, 1982, 1983;  de  Haas  et  al.  1983;  Kunst  et 
al.  1983).  Conducting  electroirs  in  these  studies 
were  excited  either  by  a  short  X-ray  pulse  of  about 
0.5  ns,  produced  by  irradiation  of  a  Pt  target  with 
3-MeV  electrons  from  a  Van  de  Graaf  accelerator, 
or  directly  by  the  3-MeV  electron  beam.  The  life¬ 
time  of  electrons  in  the  cmduction  band  is  ex¬ 
tremely  small;  at  T  =  -20®C  it  is  several  nanosec¬ 
onds  and  at  T  = -60°C  several  tens  of  nanoseconds 
(de  Haas  et  al.  1983).  The  main  processes  limiting 
the  lifetime  of  electrons  at  these  temperatures  are 
trapping  and  solvation  (hydration).  In  spite  of 
small  lifetimes,  the  analysis  of  transient  electron 
currents  in  strong  electric  fields  showed  their 
mobility  to  be  Pe  =  (25  ±  5)  x  l(Hm2  y-i  s“i  (Ver¬ 
beme  et  aL  1978,  Warman  et  al.  1980, 1982).  This 
value,  as  the  reader  must  have  noticed  before,  ex¬ 
ceeds  the  mobility  of  most  mobile  protonic  charge 
carriers  in  ice  by  more  than  four  orders  of  magni¬ 
tude  (see  Table  5).  The  mobility  of  electrons  in  ice 
is  almost  independent  of  temperature  from  -60 
down  to  -120®C. 

All  conducting  electrons  move  very  rapidly 
into  the  so-called  solvated  state.  In  the  literature 
these  electrons  are  often  referred  to  as  "solvated" 
or  "hydrated"  electrons;  these  terms  are  equiva¬ 
lent.  Such  electroiis  manifest  themselves  through 
a  specific  absorption  spectrum  with  a  maximum  at 
approximately 6800 A  and  through  a  distinct  Elec¬ 
tron  Spin  Resonance  (ESR)  line,  ^tially,  solvated 
electrons  were  discovered  in  water  and  later  in  ice. 
A  great  munber  of  publications  are  dedicated  to 
solvated  electrons  and  yet  many  of  their  proper¬ 
ties  remain  unclear.  Below,  we  shall  address  ques¬ 
tions  about  the  structure,  lifetime,  optical  proper¬ 
ties,  recombination  and  mobility  of  solvated  elec¬ 
trons  in  ice. 

The  major  reason  for  electrons  first  to  become 
localized  and  then  to  move  into  the  solvated  state 
is  considered  to  be  the  polarization  of  the  sur¬ 
rounding  medium  by  the  electron.  Owing  to  a 
large  dipole  moment,  water  molecviles  tend  to 
orient  towards  the  electron.  The  electrostatic  en¬ 
ergy  of  the  state  in  which  the  electron  is  "coated" 
by  water  molecules  oriented  towards  it  would  be 
obviously  less  than  the  initial  nonpolarized  state. 
In  water,  where  each  molecule  can  be  reoriented 


extremely  fast,  such  polarization  would  seem  to 
occur  in  a  different  fa^on  than  in  ice  in  which  the 
mutual  molecular  mentation  follows  the  ice  rules. 
Nevertheless,  dte  optical  absmrptian  spectrum  of 
the  final,  fully  relax^  state  of  the  trapp^  electnxi 
in  ice  resembles  very  closely  But  of  the  solvated 
electron  in  water,  with  only  a  small  shift  in  the 
position  of  the  absorpticm  maximum  at  1.8  eV 
being  found  upon  melting  (Nilsson  et  aL  1972, 
Shubin  et  aL  1966,  Taub  and  Eiben  1968). 

Several  of  such  spectra  in  ice  and  one  in  water 
are  shown  in  Figure  16.  The  observation  of  similar 
absorption  spectra  in  water<ontaining  media  cm 
be  consider^  as  evidence  for  the  presence  of 
solvated  electrons.  This  similarity  of  the  absorp¬ 
tion  spectra  of  solvated  electrons  in  water  and  ice 
strcm^y  suggests  the  samebasicmkroscopicstruc- 
ture  for  the  solvated  state  in  both  media.  Present 
opinicHi,  based  cm  theoretical  ccmsiderations  and 
ESR  studies  of  electrons  in  low-temperature  aque¬ 
ous  glasses  (Kevan  1981),  tends  to  favor  a  cavity 
mcxlel  for  the  solvated  electrcm,  Le.,  a  mcxlel  in 
which  the  electrcm  cxxupies  a  void  surrounded  by 
oriented  water  molecul^  This  being  the  case,  the 
initial  localizaticm  process  in  ice  must  fiierefore 
involve  trapping  at  a  vacaiuy  in  the  lattice. 

Figure  17,  adopted  from  the  paper  by  deHaaset 
al.  (1983),  shows  schematically  one  of  the  possible 
struciures  of  the  solvated  electron  in  ice.  Such  a 
defect  is  a  combinaticm  of  a  water  molecule  va- 
carKy,  a  D-defect  and  an  electron  (not  shown).  In 
this  type  of  a  defect,  ttie  electrostatic  energy  is 
lowered  by  attraction  between  the  electrcm  and 
three  protons.  A  migration  of  this  defect  must  be 
accompanied  by  ccmperative  drift  of  the  vacancy 
and  the  D-defect.  Their  mobilities  could  have  pro¬ 
vided  some  estimation  of  the  solvated  elecircm's 
mobility  in  ice.  Uidortunately,  we  could  not  find  in 
the  literature  any  reliable  data  cm  diffusion  coeffi¬ 
cients  of  either  vacancies  or  D-defecrts  or  solvated 
electrons  in  ice.  The  mobility  of  solvated  electrons 
in  water  was  measured  to  be  1.84x10-7  in2  V“i  s-i 
(Schmidt  and  Buck  1966).  This  value  can  be  ccm- 
sidered  as  an  upper  limit  for  mobility  of  solvated 
electrons  in  ice. 

Both  in  water  and  ice,  solvated  electrons  are 
annihilated  by  the  recombination  reaction  (Muto 
etal.  1992) 

e-  +  H30+^H20  +  H'  (74) 

In  addition,  they  can  react  with  impurities  in 
the  bulk  and  at  the  surface.  All  these  factors  reducre 
their  lifetimes,  which  ncmetheless  can  be  very  large 
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at  low  temperatures  (about  80  minutes  at  4.2  K 
(Muto  et  al.  1992]). 

In  eq  74  we  have  used  the  symbol "  * "  to  denote 
a  "radical,"  i.e.,  an  atom  or  a  molecule  that  has  an 
unpaired  electron  spin.  Owing  to  the  presence  of 
such  an  electron,  the  radicals  are  extremely  chemi¬ 
cally  active  and  corxsequently  unstable.  In  eq  74  we 
have  a  hydrogen  radical  or  simply  a  hydrogen 
atom  H*;  H*  is  not  the  only  radical  produced 
during  iotxization  of  water  molecules  in  ice.  We  list 
below  some  most  frequently  encountered  reac¬ 
tions  in  which  radicals  consisting  of  hydrogen  and 
oxygen  atoms  are  created  and  aimihilated: 


H2O  +  hv  — »  H20''  +  e“ 

(75) 

H2O+  +  H2O  ^  H3O+  +  OH* 

(76) 

H3O+  -i-  e-  -»  H2O  +  H* 

(77) 

OH*  +  OH*  -)  H2O2 

(78) 

H*  +  H*  ^  H2 

(79) 

H*  +  OH*  -»  H2O 

(80) 

OH*  +  H20si5H30+  +  0- 

r*4 

00 

H202  +  H*-»H20  +  0H* 

(82) 

H2O2  +  OH*  -4  HO2*  +  H2O . 

(83) 

The  study  of  the  properties  of  radicals  belongs  to 
the  extensive  domain  of  radiation  chemistry  and 
falls  outside  the  scope  of  this  report.  The  reader 
interested  in  the  properties  of  radicals  in  ice  can 
find  a  detailed  list  of  references  on  this  question  in 
the  papers  by  Bednarek  and  Flonka  (1987a,b)  and 
also  in  the  already  cited  references  concerning  free 
(conduction)  and  solvated  electrons. 

Let  us  indicate  here  that,  because  of  the  frac¬ 
tional  charge  of  ions  in  ice  (±  0.62e),  some  of  the 
reactions  (eq  75-83)  (written  for  water,  where  ions 
have  an  integer  charge)  must  be  altered  and 
completed  to  fulfill  the  charge  conservation  law. 
For  example,  the  charge  of  particles  in  the  left- 
hand  side  of  eq  77  is  -0.38e,  and  that  of  the  right- 
hand  side  is  zero.  The  charge  conservation  law  will 
hold  if  we  add  to  the  left-hand  side  a  D-defect.  In 
this  case  the  electron  is  a  solvated  electron  trapped 
by  the  complex:  a  water  molecule  vacancy  +  D- 
defect.  Thus,  now  the  reaction  (eq  74)  can  be 
written  as 

H3O  +  (e-  +  VH20  +  D)  H2O  +  H*  +  Vh20-  (84) 


For  the  charge  to  be  conserved  in  the  above  reac¬ 
tion  (eq  77),  we  have  either  to  add  a  D-defect  to  the 
right-hand  side  or  to  assume  that  the  radical  OH* 
has  a  charge  of  +0.38e.  The  possibility  for  the 
radicals  in  ice  to  be  electrically  charged  was  first 
suggested  by  Chesnakov  et  al.  (1987).  The  reasons 
for  the  radicals  embedded  in  the  ice  structure  to 
have  an  electric  charge  are  the  same  that  lead  to  the 
fractional  charges  of  ions  in  ice.  There  is  no  experi¬ 
mental  evidence  available  that  supports  or  refutes 
the  presence  of  electrical  charge  at  radicals  in  ice. 
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